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Abstract.  In  this  paper  we  prove  the  existence  and  uniqueness  of  path-wise  strong  solution  to  stochastic  viscous  flow  in 
unbounded  channels  with  multiple  outlets  using  local  monotonicity  arguments.  We  devise  a  construction  for  solvability 
using  a  stochastic  basic  vector  held. 
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1.  Introduction 

This  paper  concerns  with  stochastic  fluid  dynamics  in  unbounded  channel  domains  with  noncompact 
boundaries  generalizing  the  deterministic  results  in  Sritharan  [62] .  Mathematical  theory  of  viscous  incom¬ 
pressible  flow  through  unbounded  channel  has  many  applications  such  as  hydraulics  in  water  resources, 
hydraulic  machinery,  oil  transport  networks,  flow  in  engines  etc.  Well-posedness  theorem  is  an  essen¬ 
tial  step  for  applications  in  optimal  control  theory  (Sritharan  [64]),  convergence  of  numerical  algorithms 
and  nonlinear  filtering  (Sritharan  [63],  Fernando  and  Sritharan  [23]).  Solvability  theory  of  generalized 
solutions  to  Navier-Stokes  equations  was  pioneered  by  Leray  [41],  Hopf  [31]  and  Ladyzhenskaya  [35,36]. 
Steady  state  flow  through  channels  of  various  kinds  has  been  studied  by  a  number  of  authors  including 
Amick  [4,5],  Amick  and  Fraenkel  [6],  Ladyzhenskaya  and  Solonnikov  [37-39].  In  [4],  Amick  discussed  the 
steady  flow  of  viscous  incompressible  fluid  in  channels  and  pipes  in  two  and  three  dimensions  which  are 
cylindrical  outside  some  compact  set.  The  paper  by  Hey  wood  [28]  highlighted  the  question  of  uniqueness 
of  the  solution  of  the  Navier-Stokes  equations  for  certain  unbounded  domains  modeling  channels,  tubes, 
or  conduits  of  some  kind  and  the  importance  of  prescribing  flux  or  the  overall  pressure  difference.  In  [6], 
Amick  and  Fraenkel  studied  steady  state  solutions  of  the  Navier-Stokes  equations  in  various  types  of 
two  dimensional  channel  domains.  In  [29],  Heywood  constructed  classical  solutions  of  the  Navier-Stokes 
equations  for  both  stationary  and  non-stationary  boundary  value  problems  in  arbitrary  three-dimensional 
domains  with  smooth  boundaries.  The  time  dependent  flow  through  the  three-dimensional  channels  with 
outlets  diverging  at  infinity  has  been  studied  by  Ladyzhenskaya  and  Solonnikov  [37]  and,  Solonnikov 
[59] .  The  paper  by  Solonnikov  [59]  presents  solvability  of  boundary  value  problems  for  the  Stokes  and 
Navier-Stokes  equations  in  noncompact  domains  with  several  oulets  to  infinity.  Babin  [7]  considered  the 
Navier-Stokes  system  in  an  unbounded  planar  channel-like  domain  and  proved  that  when  the  external 
force  decays  at  infinity,  the  semigroup  generated  by  the  system  has  a  global  attractor  and  its  Hausdroff 
dimension  is  finite  using  weighted  Sobolev  estimates.  The  paper  by  Sritharan  [62]  addressed  the  following 
two  important  cases  which  were  not  considered  in  the  earlier  works: 

(i)  time-dependent  flow  through  two  and  three  dimensional  channels  with  finite  cross  section; 

(ii)  time-dependent  flow  through  two-dimensional  channels  with  outlets  diverge  at  infinity 
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Fig.  1.  Multi-channel  domain 


and  provided  a  unique  solvability  theorem  for  the  two-dimensional  case  of  the  problem  type  (i).  Solvability 
of  stochastic  Navier-Stokes  equations  in  unbounded  channel-like  domains  with  non-zero  flux  condition 
have  remained  as  an  open  problem  in  both  two  and  three-dimensions.  To  the  best  of  the  authors  knowl¬ 
edge,  this  work  appears  to  be  the  first  systematic  treatment  of  stochastic  two-dimensional  Navier-Stokes 
equations  in  such  domains.  In  this  paper  we  consider  a  stochastic  version  of  the  problem  of  type  (i)  in 
multi-channel  domains  in  two  dimensions  and  prove  a  unique  solvability  theorem  with  a  possible  future 
extension  to  three  dimensions  (up  to  a  stopping  time  determined  by  the  size  of  the  flux  and  the  Reynolds 
number).  The  problem  of  type  (ii)  may  possibly  be  resolved  by  suitably  choosing  a  conformal  mapping 
(see  Amick  and  Fraenkel  [6]  for  similar  ideas  in  the  case  of  steady  flows)  to  straighten  the  diverging 
outlets. 

Let  us  consider  the  unbounded  multi-channel  domains,  with  several  outlets  as  shown  in  the  figure  (see 
Fig.  1).  Let  the  outlets  of  the  multi-channel  domain  be  named  as  0i,02,...,0Ar  and  outside  a  compact 
region  let  the  outlets  be  of  constant  widths  di,  (i2, . . . ,  Our  first  step  is  to  construct  a  basic  vector 
field  through  each  of  these  outlets  with  the  stochastic  flux  such  that  Ylf=i  =  0?  s. 

The  methodology  of  proof  can  be  understood  by  considering  a  channel  with  two  outlets  having  a  unit 
width  connected  in  a  smooth  way  0  =  OiUO2U0o  (see  Fig.  2). 

Let  us  now  discuss  the  problems  of  type  (i)  and  examine  the  difficulties  that  arise  in  proving  solvability. 
The  time-dependent  Navier-Stokes  problem  is  usually  treated  using  the  method  of  Hopf  [31]  in  the 
deterministic  setting  which  relies  on  L^-energy  estimates.  The  traditional  methods  of  solvability  fail  in 
the  absence  of  an  energy  inequality.  For  channels  of  finite  cross  section,  as  pointed  out  in  Sritharan  [62], 
in  order  for  the  net  flux  to  be  nontrivial,  the  velocity  field  should  not  decay  to  zero  at  infinity  (upstream 
and  downstream)  and  hence  such  velocity  fields  would  then  have  infinite  energy. 


Fig.  2.  Channel  with  two  outlets  having  unit  width 
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Below  we  give  a  heuristic  argument  regarding  the  infinite  energy  of  such  velocity  fields.  We  also  point 
out  that  in  the  absence  of  a  rigorous  decay  theory,  only  a  heuristic  argument  could  be  made  in  this  regard. 

Let  the  velocity  field  be  stochastic  and  modeled  on  a  complete  probability  space  P)-  If  the 

stochastic  net  fiux  of  the  velocity  field  u  is  a;  G  then  across  any  cross  section  L,  we  have 

J  u  ■  ndS  = 

where  n  is  the  normal  to  the  curve  L  and  d^S  is  the  length  element. 

In  this  case,  the  velocity  |u(x,  ^  0  as  \x\  oo,  (where  x  =  (x^y)  with  y  is  of  constant  width) 

P-a.  s.  To  see  this  let  us  take  the  2-D  case  with  the  outlet  O2  =  {{x^y)  G  (0,  00)  x  (0,1)}.  The  fiux 
across  any  cross  section  is  same  throughout  the  channel,  due  to  divergence  free  condition.  That  is,  if 
|u(x,  t,C(;)|  ^  0  as  X  ^  00,  then  the  fiux  at  the  far  field  is  zero.  Hence  the  fiux  across  any  cross  section 
is  zero  throughout  the  channel  giving  the  net  fiux  is  zero.  Thus  for  the  fiux  to  be  non-zero,  we  need  the 
condition  that  |u(x,  t,C(;)|  ^  0  as  x  ^  00. 

As  an  example^  we  consider  the  decay  of  \x{x^y^t^uo)  in  the  following  form 

- - VC2{y,t,u}),a  >  1 

\X\ 

for  sufficiently  large  M  with  \x\  >  M,  where  Ci{y^t,uj)  and  C2{y^t^uo)  >  0  be  L^(f2;  C(0,  T;  L^(0, 1))) 
functions  such  that 

/  Ci{y,t,Lo)dy  =  Q  and  f  C2{y,t,u})  dy  = 

Jo  Jo 

From  above,  we  have  u(x,  y,  t,  lo)  is  fiux  carrying  in  the  far  field  as 

[  u(a:,y,t,w)dy  ~  dy  +  f  C2{y,t,u))dy  = 

Jo  Jo  1^1  Jo 

Then  at  any  given  time  t  G  (0,T],  the  L^-energy  is  given  by 

nOO  pOO 

\u{x,y,t,uj)\'^  dxdy  >  /  /  \u{x,y,t,uj)\'^  dx  dy 

Jo  Jm 

Ci{y,t,oj) 


/o  Jm 

nl  pOO 


+  C2{y^  L  ^) 


JO  JM 

pi 

+  2 

Jo  Jm 


kr“  Jo  Jm 

Ci{y,t,u>)C2{y,t,u>) 


dx  dy 

\C2iy,t,u))f  dx  dy 


dx  dy. 


(1.1) 


The  first  integral  in  (1.1)  always  has  a  finite  positive  value  and  is  equal  to 

(2o  -  1)AP°-'  (i 

The  second  integral  diverges  to  Too.  The  last  integral  in  (1.1)  is  also  bounded,  since 
Ci{y,t,oj)C2{y,t,uj) 


aoo 

4 


X\ 

1  pOO 


< 


< 


dx  dy 


\Ci{y,t,uj)\\C2{y,t,uj)\ 


dx  dy 


0  Jm 


{a  -  1)M' 


)M«-1  (/  \C2{y,t,uj)\^dy 


1/2 
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Thus  for  any  given  time  t  >  0,  we  have 

oo 

|u(x,^,  =  oo,  P-a.  s. 

There  are  extensive  literature  on  deterministic  flow  through  channel  type  domains.  Interested  readers 
may  look  into  Amick  [4,5],  Amick  and  Fraenkel  [6],  Babin  [7,8],  Borchers  and  Pileckas  [12],  Heywood 
[28-30],  Kapitanskh  and  Piletskas  [32],  Ladyzhenskaya  and  Solonnikov  [37-39],  Pileckas  [52],  Piletskas 
[53],  Solonnikov  [59],  Solonnikov  and  Piletskas  [60],  Sritharan  [61,62]  to  name  a  few. 

For  a  sample  of  literature  on  stochastic  Navier-Stokes  equations,  we  refer  the  readers  to  Bensoussan 
[9],  Bensoussan  and  Temam  [11],  Capinski  and  Cutland  [14],  Da  Parto  and  Zabczyk  [18],  Flandoli  and 
Gatarek  [24],  Menaldi  and  Sritharan  [47],  Pardoux  [50,51],  Sritharan  and  Sundar  [65],  Vishik  and  Fursikov 
[68],  Sritharan  [64],  Fernando  and  Sritharan  [23],  Sakthivel  and  Sritharan  [55]. 

The  plan  of  the  paper  is  as  follows.  In  Sect.  2,  the  main  result  of  this  paper  and  the  functional  setting 
have  been  given.  A  divergence  free  vector  field  of  infinite  energy  carrying  a  nontrivial  net  flux  through 
the  channel  is  constructed  in  Sect.  3  using  the  solution  of  the  heat  equation.  In  Sect.  4,  we  characterize 
the  properties  of  the  linear  and  bilinear  operators  that  are  associated  with  the  Navier-Stokes  problem.  A 
perturbed  vector  field  is  constructed  in  Sect.  5  using  a  suitable  transformation  involving  the  constructed 
basic  vector  field.  A-priori  estimates  for  the  solutions  of  the  perturbed  vector  field  are  obtained  in  Sect.  6. 
In  Sect.  7,  we  prove  the  local  monotonicity  condition  for  the  sum  of  the  Stokes  and  the  inertia  operators 
as  well  as  the  existence  and  uniqueness  of  strong  solutions  to  the  perturbed  vector  field  by  exploiting  this 
local  monotonicity  condition.  In  Sect.  8,  we  mathematically  characterize  the  perturbation  pressure  field 
using  a  generalization  of  the  de  Rham’s  Theorem  to  processes.  Section  9  completes  the  proof  of  the  main 
result. 


2.  Basic  Definitions  and  the  Main  Theorem 

In  this  section,  following  Sritharan  [62] ,  we  define  the  class  of  channel  domains  that  will  be  analyzed. 

Definition  2.1.  (Admissible  channel  domain)  A  simply  connected  open  set  0  C  with  boundary 
90  consisting  of  two  disconnected  components  90 1  and  902  is  called  an  admissible  channel  domain  (see 
Fig.  3),  if  it  is  the  union  of  three  disjoint  sets  0o  U  Oi  U  O2  defined  in  the  foiiowing  way.  Let  Oi  and 
O2  be  two  semi- infinite  strips  of  width  di  and  d2  respectively.  These  two  straight  channels  are  smoothly 
(not  necessarily  coaxially)  joined  by  a  bounded  domain  0o  such  that  90i  U  902  =  90  G  C^. 

Now  let  us  consider  the  problem  of  accelerating  a  viscous  incompressible  fluid  from  rest  to  a  given 
stochastic  flux  rate  through  an  admissible  channel  domain.  Let  (f2, .F,  .Ft,  P)  be  a  complete  probability 
space.  The  mathematical  problem  is  to  find  the  velocity  field  u  and  pressure  field  p  such  that 

(u,p)  :  0  X  [0,  T]  X  ^  X  M, 

the  momentum  equation 

Ut  +  u  •  Vu  =  —  Vp  +  uAu  +  ^(x,  t,  cj)  in  0  X  (0,  T)  x  (2.2) 


Fig.  3.  Admissible  channel  domain 
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the  incompressibility  condition 


V  •  u  =  0  in  0  X  (0,  T)  X  12, 

(2.3) 

the  non-slip  boundary  condition 

on  the  channel  walls 

the  initial  condition 

u(x,t,u;)  =  0  on  30  x  [0,T]  x  12, 

(2.4) 

u(x,  0,^;)  =  0,  (x,  cj)  G  0  X  12, 

(2.5) 

and  the  flux  condition 

J  u  •  nd5'  =  for  all  t  G  [0,T]  with  ^(0,0;)  =  0,  V  cj  G  (2.6) 

are  satisfied.  The  properties  of  the  stochastic  flux  will  be  discussed  in  the  later  sections.  Here  >  0  is 

the  coefficient  of  kinematic  viscosity  and  T  is  any  cross-sectional  curve  cutting  the  channel. 

In  this  formulation  the  stochasticity  of  fluid  flow  is  due  to  an  external  random  forcing  and  the  random 

flux.  Also  we  will  assume  that  the  external  random  forcing  ^(x,  t,^;)  and  the  random  flux  are 

mutually  independent  processes.  Further  details  about  the  noise  have  been  provided  in  the  subsequent 
sections. 

Now  let  us  state  the  main  result  of  this  paper. 


Theorem  2.2.  Suppose  that  the  flux  rate  satisfies  the  moment  bound: 


<  Ci(r), 


(2.7) 


for  some  preseribed  T.  Then  for  eaeh  sueh  there  exits  a  unique  strong  solution  u(x,  t,a;)  with  the 

following  estimates: 


E 


[  sup  [  |(u  —  w)(x,  t,  •)|^  dx  +  2z/  /  [  |V(u  —  w)(x,  t,  •)|^  dxdt  ) 

\o<t<T  Jq  Jo  Jq  J 


(2.8) 


for  some  divergenee  free  veetor  field  w(x,  t,  uo)  that  vanishes  on  the  boundary  30  and  earries  the  preseribed 
flux 


■  ndS  =  ^{t,uj),  V  a;  G  12. 


(2.9) 


We  prove  the  above  theorem  in  the  subsequent  sections.  The  following  functional  frame  work  is  used 
in  this  paper. 


(7^(0)  =  The  space  of  all  infinitely  differentiable  vector  fields  with 
compact  support  in  0, 

1^0  (0)  =  The  completion  of  (7^(0)  vector  fields  in  the  norm  ||  V(/:?||l2(©), 

JJ1(0)  =  {<^  ;  0  ^  m2.  ^  g  L2(0)}, 

Hl{@)  =  {^  :  0  ^  R2;  ^  G  L2(0),  G  L2(0)  and  =  O}, 

j(0)  =  :  0  ^  m2;  ip  G  C^(0)  and  V  ■  <p  =  O}. 
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Also  let  US  define 

H  =  Completion  of  j(0)  in  the  L^(0)  norm, 

V  =  Completion  of  j(0)  in  the  H^{Q)  norm, 

H*  =  :  0  ^  cp  G  L^(0);  V  •  (/?  =  0  and  cp  •  n|^Q  =  O}, 

V*  =  {ip:e^R^;  ipe  and  V  •  (/^  =  O}, 

Vo  =  Completion  of  j(0)  in  the  norm  of  ||  V(/?||l2(©), 

Vo*  =  {(/^GlVo(0);V.^  =  O}. 

Let  us  denote  the  norm  in  H  by  |  •  |  and  the  norm  in  V  by  ||  •  ||.  If  we  identify  H  with  its  dual 
H'  using  the  Riesz  representation  theorem,  we  obtain  the  continuous  and  dense  embedding 

V  C  H  =  H'  C  V'. 

Also  let  us  denote  the  duality  pairing  between  V  and  V'  by  (•,  •).  Note,  however,  that  (unlike  in  bounded 
domains)  the  embedding  V  C  H  is  not  compact  since  0  is  unbounded.  Poincare  lemma  holds  for  admis¬ 
sible  channel  domains  (since  they  have  finite  cross  section): 

||0||L=(e)  <C||V0||L=(e),  V,/.Gffo'(©) 

and  hence,  in  V  the  norm  of  is  equivalent  to  that  obtained  by  the  Dirichlet  integral  ||  V(/)||l2(©). 


3.  Construction  of  the  Basic  Vector  Field 

Following  the  ideas  from  Sritharan  [62],  we  will  construct  a  divergence  free  basic  vector  field  w(x,  t,^;) 
in  0  which  vanishes  on  50  and  carries  the  prescribed  random  flux  through  the  channel.  Note 

that  this  “constructed”  vector  field  need  not  satisfy  the  Navier-Stokes  equations  in  0  although  it  does 
in  Oi  and  O2  due  to  the  nature  of  the  construction  used.  The  method  can  be  described  as  follows.  Using 
the  solution  of  the  one-dimensional  heat  equation  with  random  flux  ^(t,a;),  a  vector  field  is  constructed 
in  each  of  the  straight  channel  outlets  Oi  and  O2.  We  then  smoothly  join  these  two  vector  fields  by 
constructing  a  smooth  extension  in  0o.  In  Oi  and  O2  the  vector  field  will  have  only  one  component, 
namely  in  the  direction  of  the  axes  of  the  outlets.  In  0o,  however,  in  general  both  components  of  w(x,  t,  uj) 
will  be  nonzero.  Let  us  now  consider  one  of  the  outlets,  say  O2  and  assume  for  simplicity  that  the  width 
of  the  channel  is  unity  (see  Fig.  4). 


y 


A 


X 


Fig.  4.  The  outlet  O2 
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Let  us  define  the  outlet  O2  =  {{x^y)  G  (0,  00)  x  (0, 1)}.  Let  us  seek  a  divergence-free  vector  field  in 
the  form  t,  a;),  0)  and  a  scalar  field  P{x,t^uj)  in  O2  such  that 

d  d 

—wi{y,t,uj)  -  u—wi{y,t,uj)  =  -  —  P{x,t,uj)  =  in  O2  x  (0,T)  x  ft,  (3.1) 

wi{0,t,Lu)  =  wi{l,t,Lu)  =0  for  t  e  [0,T],cj  G  (3.2) 

u;i(^,  0,  a;)  =  0  for  ^  G  (0, 1),  a;  G  and  (3.3) 

P 

Wi  (^,  L  uo)  dy  =  ^(t,  u)  for  t  G  [0,  T]  with  ^(0,  cj)  =  0  V  cj  G  fl.  (3.4) 


dt 


/ 


Here  the  function  /(t,  u)  needs  to  be  determined  from  the  prescribed  flux  ^(t,  uo).  To  resolve  this  problem 
we  first  write  down  the  solution  of  the  system  (3.1)-(3.3)  in  terms  of  f{t^uo)  and  then  use  the  condition 
(3.4)  to  evaluate  f{t,uj)  in  terms  of  The  solution  of  (3.1)-(3.3)  can  be  obtained  by  the  method 

of  separation  of  variables  (see  Cannon  [17]).  The  existence  and  uniqueness  of  solution  of  the  boundary 
value  problem  for  the  heat  equation  with  stochastic  boundary  conditions  has  been  proved  in  Cahlon  [16]. 

Theorem  3.1.  Let  f{t^uj)  satisfies  the  following  uniform  Holder  rendition  on  [0,T]^  for  all  tiH2  ^  [0^7^]; 


<  L{uj)\ti  -  t2r,o  <  7  <  2’ 


(3.5) 


with  L{uj)  >  0  and  L(-)  G  L^(f2).  Then,  there  exists  a  unique  solution  wi{-,  ■,  •)  of  the  problem  (3.1)-(3.3) 
sueh  that 

wi{y,t,io)eL^  (^^;C([0,1]  X  [0,r])) 

and  satisfying  the  following  eonditions. 

(i)  wi{yH^^)  is  a  eontinuous  funetion  of  y  and  t  in  L^(fl).  i.e., 


{wi{y,t  +  h,uj)  -Wi{y,t,Lo)y 
{wi{y  +  h,t,uj)  -  wi{y,t,u>)y 


0  as  h  ^  0, 
0  as  h  ^  0. 


(ii)  There  exists  a  stoehastie  funetion  wi^{yHi^)  ^  L^(^;  ^^([0?  1]  x  [0^7^]))  sueh  that 


E 


Wi  {y,  t  +  h,u))-wi  {y,  t,  oj) 


-  u>u(h,t,aj) 


0  as  h  ^  0. 


(hi)  There  exists  a  stoehastie  funetion  wy  (^,  t,a;)  G  L^(fl;  (7([0, 1]  x  [0,T]))  sueh  that 


E 


wi{y  P  Pwi{y  -  h.t^uj)  -  2u;i(^,t,cj) 
h^ 


0  as  h  0. 


(3.6) 


(3.7) 


(3.8) 


(iv)  The  equation  wi^{y Hi  =  ^^lyyivPiX))  +  /(t,c<;)  is  satisfied  for  a.e.  a;  G  Ll. 

Proof.  By  using  the  method  of  separation  of  variables,  we  can  write  down  the  solution  of  (3.1)-(3.3)  as 

wi{y,t,u})  =  -'p{f  sin(2n  +  l)7ry.  (3.9) 

2n  +  l  J 


We  now  use  the  Weierstrass’  M-test  to  show  that  the  infinite  series  solution  (3.9)  to  the  Heat  Equation 
(3.1)-(3.3)  converges  uniformly.  Let  us  take 

an{y,t,uj)  =  f  /(s,  w)— 2— +  l)7r2/,  for  t  >  0. 

Jo  2n  +  1 
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Then  by  using  Holder’s  inequality,  for  t  G  (0,T]  and  y  G  [0, 1],  we  have 

-'^(2"+!)'’^'*  |/(s,w)|2ds^  e2'^(2"+l)'^'Ms^ 


< 


(2n  +  l) 


1/2 


<  ll/(•,w)l|L2(0,T) 


^g2zy(2n+l)^7r^t  _  ^ 

2u{2n  +  1)^71^ 


1/2 


(2n  +  l)^7r\/^^ 


-a.  s. 


(3.10) 


Hence,  we  get 


|an(2/,i,w)|  <  M„  where  M„  =  ||/(•,  w) ||l2(o,t) - 


(2n  +  l)^Tr-\/^’ 

Since  X;“=o  (2^^4tf  =  ^  ^rom  Apostol  [2],  we  have  X;“=i  then  X:“=i  +  E“=o  (SiltlF  = 

^  gives  Er=o  (2n+i)"  =  J2n=o^n  =  ^^\\f{-,oj)\\L^(o,T),  P-a.  s.,  for  all  t  e  (0,T]. 

Hence  by  the  Weierstrass’  M-test,  wi{y,t,u>)  =  ^J2^=QO.n{y,t,u>)  converges  absolutely  and  uniformly 
for  y  G  [0, 1]  and  t  >  0.  For  t  =  0  also  the  convergence  is  uniform  since  wi{y^  0,a;)  =  0. 

One  can  re-write  (3.9)  as 


where 


From  (3.5),  we  have 


wi{y,t,uj)=  [  K{y,t-s)f{s,uj)ds={K*f){t), 
Jo 

.  oo  ^ 

K{y,  t)  =  -J2  sin(2n  +  l)7ry. 

^G;,  (2n  +  i) 


(3.11) 


E  [|/(ti,u;)  -  fit2,uj)f]  <  E  [L(u;)|ti  -  ^211' 

=  E\L{io)\^\ti  -  ^  0  as  ti  ^  t2, 


for  all  ti,t2  G  [0,T],  which  gives  /(*,•)  G  L^(f2;  C[0,  T]).  We  can  prove  that  the  solution  (3.9)  of  the 
problem  (3.1)-(3.3)  satisfies  wi{y,t,uj)  G  (^2;  C([0, 1]  x  [0,T]))  (see  Cahlon  [16]  and  Walsh  [69]  for 
existence  and  uniqueness). 

Let  us  now  prove  part  (iv)  of  the  theorem.  For  t  >  0,  we  have 


d 

Wt 


{wi{y,t,uj)) 


[  K{y,t-  s)f{s,uj)ds 
Jo 

nOO 

/  H{t  -  s)K{y,t  -  s)f{s,u})ds 
Jo 


dt 
dt 

f 

f)  A  ^ 

/  s)f{s,u>)ds  +  -J2  ' 

Jo  dt  tt 

J  ^-^(2/,t-s)/(s,w)ds  +  /(t,w), 


H{t  -  s)—K{y,t-  s)  +  K{y,t-  s)5{t  -  s) 


sin(2n  +  l)7ry 


(2n  +  1) 


/(s,w)ds 


(3.12) 
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where  H{')  is  the  Heaviside  function.  The  last  term  is  obtained  by  using  the  property  of  (^-function  and 
using  the  formula  ^^^(2n+i)^^  ~  f’  ^  ^  By  a  direct  calculation,  as  is  uniformly 

convergent  in  time  t  >  0  for  a  fixed  y  G  (0, 1),  it  is  easy  to  see  that 

^K{y,  t-s)  =  -4z/7r  y^(2n  +  i)e-^(2"+i)"^"(t-«)  sin(2n  +  l)7ry. 

n=0 


From  (3.12),  using  the  uniform  convergence  of  the  above  series  in  time  t  >  0  for  a  fixed  y  G  (0, 1),  we 
have 


d 


—  (M)i(2/,i,w))  = -4j/7r^  <{  I  (2n  +  l)e  '^(2"+i)  ^  «)/(s,w) 


r*  r 

n=0  '"^0 

X  sin(2n  +  l)7ry  |  +  f{t,  co). 

Next,  let  us  calculate  -^wi{y,t,u)).  We  have 


ds 


(3.13) 


/■* 

-^wi{y,t,u})  =  -^[K{y,t  -  s)]f{s,u))ds.  (3.14) 

Since  K{-,-)  is  uniformly  convergent  in  y  and  ^(-j  •)  is  also  uniformly  convergent  in  y  for  a  fixed  t  >  0, 
by  a  simple  calculation,  we  have 


9y2 


[K{y,t 


-s)] 


-47r  y](2n  +  i)e-‘^(2"+i)"^"(*-«)  sin(2n  +  l)7ry. 

n=0 


The  above  series  is  uniformly  convergent  in  time  t  >  0  for  a  fixed  y  €  (0, 1)  and  hence  from  (3.14),  we 
get 

f\2  ^  rt 

=  — 47r  /  (2n  +  ^  a;)  sin(2n  +  l)7r^d5.  (3.15) 

t'o-Jo 

The  above  integral  is  well  defined  and  from  Eqs.  (3.13)  and  (3.15),  we  have  part  (iv)  of  the  theorem. 

Now  we  prove  part  (i)  of  the  theorem.  Since  K  is  continuous  in  time  t,  for  any  given  6  >  0,  there 
exists  a  77  >  0  such  that 


\K{y,ti)  -  K{y,t2)\  <  s  whenever  |ti  -  ^2!  <  ri, 

ti,t2  G  (0,  T].  For  /(•,•)  G  L^(f2;  C[0,  T]),  by  choosing  \h\  <  7;,  let  us  consider  E[u;i(7/,  t+h,  Lo)—Wi{y,  t,  lo)\^ 
and  by  using  Young’s  inequality  and  continuity  of  K {y,  t)  in  t,  we  have 

E  [wi{y,t  +  h,u!)  -  wi{y,t,u>)f 


=  E^y  K{y,t  + h- s)f{s,Lo)ds  -  J  K{y, t  -  s)f{s, u>)d^ 

<  2E  \K{y,t  +  h  -  s)  -  K{y,t-  s)|^  |/(s,w)pds^ 

+  2E  \Kiy,t  +  h  -  s)f\fis,L0)\^  dsj 

<  2E  sup  |/(5,c<;)p  f  \K{y,t h  —  s)  —  K{y,t  —  s)\‘^  ds 

_0<s<t  _  Jo 
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+  2E 


sup  |/(s,w)| 

P<S<t 


h  —  s)\‘^  ds 


<mf 


L2(O;C[0,T]) 


/t+n 

\Kiy,i 

^  \K{y,t  ^  h  —  s)\‘^  ^  0  as  ^  0, 


since  5  is  arbitrary  and  as  ^  0,  t  ^  5  and  K{y,t  h  —  s)  ^  K{y,0)  =  1  for  all  y  G  (0, 1).  Similarly, 
'E{wi{y  +  /i,  t,  Lu)  —  wi{y^  t,  cj))^  — ^  0  as  ^  0. 

Let  us  now  prove  part  (ii)  of  the  theorem.  Since  for  a  fixed  y  G  (0, 1),  K{y,t)  is  uniformly  convergent 
in  time  t  >  0  and  its  derivative  ^  exists  and  also  is  uniformly  convergent  for  t  >  0,  we  have  for  a  given 
5  >  0,  there  exists  an  77  >  0  such  that 


K{y,t  +  h)-K{y,t) 
h 


<  e  when  ever  \h\  <77. 


To  prove  (3.7),  let  us  use  the  differentiability  of  K{y^t)  in  time  t.  For  \h\  <  77,  we  have 


E 


wi{y,t^  h,u)  -  wi{y,t,Lo) 


h 


K{y,t  +  h-  s)  -  K{y,t-  s) 
h 


f{s,uj)ds 


a 

1  /■*  d 

^hj  '^(2/’*  + ^  -  7  -^K{y,t  -  s)f{s,u})ds  -  f{t,u}) 


<  2E 


sup  \f{s,u))\ 

,0<s<t 


1/ 


K{y,t  +  h- s)  -  K{y,t- s)  d 
- h - 


ds 


+  2E  - 


< 


/  ^  nt+h 

{-J. 

‘^\\f\\hi(i-,c[o,T]f^T  +  2E  (  t  ^  [K{y,t  +  h-  s)/(s,  w)  -  ds 


,t  +  h-  s)f{s, u>)  -  f(t, w)]  ds 

pt-\-h 


Note  that  in  the  time  interval  [t,  t  +  h]  as  h  ^  0,  s  ^  t  and  hence  the  function  K{y^t  ^  h  —  s)f{s,  uo) 
K{y,0)f{t,uj)  =  Thus  by  the  Lebesgue’s  differentiation  theorem  (Theorem  6,  Appendix  E.4  of 

Evans  [21]),  the  last  term  of  the  right  hand  side  of  the  above  inequality  goes  to  0  as  h  ^  0.  Einally,  since 
5  >  0  is  arbitrary,  we  have  the  desired  result  (3.7).  Similarly  one  can  prove  that  there  exists  a  stochastic 
function  (^,  t,  cj)  G  L^(fl;  C([0, 1]  x  [0,T]))  such  that 


E 


wi{y  +  h,t,uo)  ^wi{y  -  h,t,uo)  -  2wi{y,t,uo) 


-wi„{y,t,uj) 


0  as  h  — >  0. 


□ 


Corollary  3.2.  For  G  L^(fl;  C[0, T]),  there  exists  a  funetion  ^^(t)  ^  sueh  that 

^^{t,u})  =  +  J  ^[)(t  -  r)/(r,w)dr, 

where  i){t)  =  ^  (2;4tF  exp(-;^(2n  +  ifTTH). 


(3.16) 
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where 


..  „„.J«  (2»  +  l)^ 

=  / 

Jo 

8  1 

=  (2n  +  ir 

n=0  ^  ' 


(3.17) 


(3.18) 


Note  that  since  \){t)  is  a  decreasing  function  of  t,  we  have  for  all  t  >0,  mt)  <m  =  =  ^- 

Also  the  series  [)(t)  is  convergent.  Then  there  exists  a  function  G  C[0,T])  such  that 


dt 


=  f{t,uj)  +  ^[l(t  -  T)/(T,w)dT, 


(3.19) 


where  =  — 8z/ e  ^(2n+i)^7r^t^  Since  ()  is  differentiable  in  time  t  >  0,  for  any  given  5  >  0,  there 

exists  an  7^  >  0  such  that 


[)(t  +  h)  -  [)(t)  d 

— h — 


<  5  when  ever  \h\  <77. 


For  proving  the  existence  of  let  us  use  /  G  C[0, T]),  the  differentiability  of  ()  and 

choose  \h\  <  77  to  get 

IT. n.  (jn  —  IT. hn  n  ^ 

E 


\^{t  +  h,L0)  -  9 

h 


<  2E 


+2E 

<  2E 

+2E 


1  rt-\-h 

-  J  H)(t  +  h-  T)f{T,  U))  -  f{t,  w)]  dr 


sup  |/(t,w)|^ 

.0<r<t 

nt-\-h 


f 


dr 


1 

-  J  H)(t  +  h-  T)f{T,  u))  -  f{t,  w)]  dr 


1  rt-\-h 

-  J  [[){tPh-  r)f{T,  uj)  -  f{t,  (j)]  dr 


Note  that  in  the  time  interval  [t,  t  +  h]  as  h  ^  0,  s  ^  t  and  hence  the  function  [)(t  +  h  — T)/(r,  cj)  ^  /(t,  cj). 
Thus  by  the  Lebesgue’s  differentiation  theorem  (Theorem  6,  Appendix  E.4  of  Evans  [21]),  the  last  term 
of  the  right  hand  side  of  the  above  inequality  goes  to  0  as  h  ^  0.  The  arbitrariness  of  5  gives  the  required 
result.  □ 


We  have 


dt 


r{t,u})  =  f{t,u))  +  ^[)(t  -  r)/(r,w)dr, 


(3.20) 


d 


-m  =  -svY.' 


—  u(2n+lp7T^t 


n=0 


where 
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From  (3.20),  denoting  —  r)  as  —  r),  we  obtain  the  Volterra  integral  equation  of  the  second  kind 
for  the  determination  of  /(•,  •)  in  terms  of 

~  d'T-  (3.21) 

Lemma  3.3.  Let  ^  L^(f2;  L^(0,  T))  be  given.  Then  there  exists  a  unique  solution  /(*,•)  G 

L^(f2;  L^(0,  T))  for  the  integral  equation  (3.21)  in  the  form 


oo 


E 

n=0 


,d^ 

dt 


(3.22) 


where  is  given  by 

[J(Lfj]{t^uj)  =  —  f  7i{t  —  T)'ip{r,uj)  dr  = —{H  ^  fj){t,uj).  (3.23) 

Jo 

Proof  Let  us  begin  with  the  convolution  (3.23)  and  apply  Young’s  inequality  for  convolutions  to  obtain 

ll‘^l^(->-)llL2(n;L2(0,T))  =  E  -  J^'H{t-T)lp{T,Uj)dT  dtj 

=  E  *  '0)(',  •)IIl2(0,T)) 

<  E  ^||(H(-)||li(o,t)  ')IIl2(o,t)) 

=  ll^(')  IIl1(0,T)  ^  (ll'^('5  ')IIl2(0,T)) 

=  ll^(’) IIli(o,t)  ’) IIl2(0;L2(o,t))  (3.24) 

Let  us  define 

P  =  ll^(•)llLqo,T)  <  ^)(0)  -  m  =  1  -  m  <  1- 

Hence  from  (3.24),  we  get 

L2(0;L2(0,T))  ^  p\mmn;LH0,T))^  VV’GL2(a;L2(0,T)).  (3.25) 

Let  us  define  the  map  ^  to  be 

\Jf\  {t,oj)  =  L^(t,uj)  +  [Jf /]  (l,w). 

Let  us  denote  the  approximate  solutions  of  (3.21)  as  /o,  /i,  •  •  • ,  /n  and  they  are  given  by 
/o(i,w)  = 


~  J  Wt-  'r)/o('r,w)dT 


d 


1  r 


—  ^{t,LO)  +  [/C/o]  it,LO)  =  i^fo]  {t,LO)  =  y] 


i=o 


dt 


/2(i,w)  =  -  J  n{t-  T)/i(T,w)dT 


(3.26) 
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d 


—  +  [^/o]  +  [/C^/o] 


dt 

= 

j=o 


dX 

dt 


fnit,U})  =  +  [ICfn-l]  {t,U})  =  [J?  fn-l] 


Clearly  these  approximations  satisfy 

fn{t,U>)  =  [ 

and 

fn{t,Uj)  -  /n-l(t,Cj)  = 

Hence  by  using  the  estimate  (3.25)  on  (3.29),  we  obtain 

||/n(b  ’)  ~  /n-l('5  ')  IIl2(O;L2(0,T))  ^  P 


fo]  it,  ^)  =  XI 

i=o 


xt 


dX 


{t,U}) 


dX 


dX 


dt 


L2(0;L2(0,T)) 


with  p  <1. 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


The  existence  of  a  unique  fixed  point  f{t,uj)  G  L^(f];  L^(0,T))  to  the  map  ^  is  ensured  by  Banach  con¬ 
traction  mapping  theorem  and  using  the  estimate  (3.30).  The  unique  fixed  point  /(t,  uj)  G  L^(f];  L^(0,  T)) 
to  the  map  ^  is  such  that 


lim  ||/(•,•)  /n('c) IIl2(0;L2(o,t))  ~  ^ 


and 


i^f]  it,^^)  =  fit,U})  =  X 


n=0 


X' 


.dX 

dt 


{t,L0). 


(3.31) 

(3.32) 

□ 


Let  us  now  go  back  to  the  system  (3.1)-(3.4)  and  note  that  a  given  ^=^(-,-)  €  L^(a;L^(0, T)) 
corresponds  uniquely  to  a  solution  tci(-,  •,  •)  given  by  (3.9)  and  (3.22).  Moreover,  P(-,  •,  •)  is  given  by 

P{x,t,oj)  =  -xf{t,oj)  +gi{t,(j), 

for  some  arbitrary  scalar  function  gi{t,u;)  (which  can  be  set  to  zero  for  a.  e.  Q)  for  each  uj  G  Cl.  The 
regularity  of  u;i(-,  •,  •)  can  be  verified  by  estimating  (3.9). 

Let  us  now  obtain  a  stochastic  version  of  the  Theorem  2  (section  3)  from  Sritharan  [62]. 

Theorem  3.4.  Let  ^^(•,  •)  G  L^(f2;  L^(0,  T))  be  given.  Then  the  solution  (wi^P)  o/ (3.1)-(3.4)  is  unique 
and  for  0  <  e  <  | 

w,  eL\n-,L‘^{o,T-,  H^io,i)  n  H^+%0,1))  nC{[o,T];H^io,m 


and 


dtWieL\fl;L\0,T;H"iO,m- 
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Proof.  We  only  need  to  verify  the  regularity.  Let  us  denote  by  ^  the  Friedrich’s  extension  of  the  operator 
—dy  to  iLo(0, 1)  n  iL^(0, 1).  Then  from  (3.9),  by  using  the  uniform  convergence  of  the  series  solution 
Fubini’s  theorem  and  Young’s  inequality  for  convolutions,  we  have 


E 


^-u{2n+lp7T^{t-T) 


L2(0,l) 


X  [sin((2n  +  l)7r^)(2n  +  l)27p2j(i+£/2) 

<  C'(7r)E  f^(2n  +  1)^+2^  y  J  /(r,  w)e-'^(2n+i)V^(t-T) 

<  C'(7r)E  (^(2n  + 1)2+2^  T /(r,  w)e-'^(2n+i)V=(t-r)  ) 

Vn=0  L2(C 

/  oo  ^ 

<C(7r)E  ^(2n  +  l)2+2^||/(.,0ll!.(o,r)l|e“'^^^”+'^'"'l^(o,r) 

\n=0 

oo  /  1 

<  C(7r)E  (||/(■,  OIIl^co.t))  ^(2^^  + 


L2(0,T)^ 


n=0 


^u{2n  +  l)^7r^ 

<  Ci(7r,T,  z/)||/(-,  •) IIl2(0;L2(o,t)) 

n=0 

Thus  noting  that  +  i'^-2-\-2s  <  oo  for  5  <  |,  we  obtain 


<  00. 


(3.33) 


(3.34) 


Then  one  can  deduce  that,  G  L^(r2;  L^(0,  T;  L^(0, 1)))  which  implies  wi  G  L^(f2;  L^(0,  T;  H^+^(0, 

1))).  Since  ^  is  the  Friedrich’s  extension  of  the  operator  —dy  to  i^o  (0, 1)  Pi  i^^(0, 1)  and  from  (3.34),  we 
get 

rT 


E 


11^2(0,1)  di 


<  OO, 


□ 


which  gives  dfWi  G  L^(f2;  L^(0,  T;  iL^(0, 1)))  for  0  <  e  < 

Let  us  now  describe  a  method  to  extend  the  constructed  flux  ^(^  •)  carried  in  into  the  domain 
©0  in  a  smooth  manner.  Let  ©o  be  an  open  subset  of  ©  with  compact  closure  such  that  ©0  C  ©0  (^  ©• 
Moreover,  9©o  is  of  class  (see  Fig.  5). 

Proposition  3.5.  There  exists  a  veetor  field  w  :  ©o  x  [0,T]  x  in  the  form  w(x,  t,^;)  = 

{fjy^  — '0a,)(x,  t,  o;)  for  some  stream  funetion  sueh  that 

V  •  w  =  0  in  ©0  X  (0,  T)  X 

in  the  sense  of  distributions  and 

w(x,  t,  cj)  =  0,  t  G  (0,  T),  X  =  (x,  y)  G  90o  H  9©, 

w{x,t,u))  =  {w['\y,t,u}),0),  t  e  (0,T),  x  =  {x,y)  e  dOo  n  Oi,  i  =  1,2, 
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for  almost  all  to  G  fl,  in  the  sense  of  trace.  Here  denotes  the  solution  of  the  heat  equation 

constructed  in  O*.  Moreover, 


w  e  L2(n;  l2(0,  T;  H^{Qo))  n  C([0,  T];  i/i(0o))) 
with  dtw  e  L2(n;L2(0,T;L2(eo)))  C  L2(n;  l2(0,  T;  iJ-i(eo))). 


Proposition  3.5  can  be  deduced  from  the  well-known  method  of  extending  the  boundary  data  (with 
zero  net  flux)  into  the  domain  as  a  divergence-free  vector  field  (Ladyzhenskaya  [36]).  Also  we  note  that 
the  regularity  of  the  boundary  data  corresponds  to  that  obtained  in  the  solution  of  the  heat  equation 
and  hence  we  have 


e  L2(n;  l2(0,  T;  Hi+^deo))  n  C([0,  T];  if  5  (SGo))) 

and 


Let  us  use  w(x, t,a;)  in  ©o  and  Wi\y,t,uj)  in  O^,  i  =  1,2  to  construct  the  desired  vector  field  w(x, Lcj) 
in  0  (see  Sritharan  [62]  for  more  details  about  the  construction). 

Let  t,  a;)  be  the  stream  function  corresponding  to  the  vector  field 


in  O^.  Let  us  set  =  0  on  the  lower  wall  of  O^,  hence  on  the  upper  wall  of  O^,  we  obtain 

f  Wi{y,t,oj)dy=  [  ^i)^^\y,t,u))dy  =  i)^^\y,t,u))  =  i!^''\l,t,u}). 

Jo  Jo  oy  0 

Also  due  to  the  flux  condition,  we  have 


^lae.naeo  =  ^lae.naOo  =  O' 

Let  us  construct  a  function  A(-)  :  M  ^  (O5 1)  by  mollifying  a  step  function  such  that 

0  for  ^  <  0, 


A(e) 


{? 


for  0  <  5o  <  ^, 


where  eo  is  a  small  number.  Let  us  define 

jlj{x,y,t,Lu)  =  (1  -  X{x))'ip{x,y,t,uj)  -h  X{x)'ip^^\y,t,uj)  in  ©o  H  0^,2  =  1,2, 


for  each  uj  E  Q. 

In  this  way  we  obtain  a  stream  function  which  takes  J^{x^y^t,uj)  in  ©0  and  smoothly  become 
t^uj)  in  Oi.  On  the  lower  wall  90i,  we  have 

'lp{x,t,U>)  =  'tl>{x,t,U!)  =  =  0 


62 


U.  Manna  et  al. 


JMFM 


and  on  the  upper  wall  502,  we  have 

Hence  we  obtain  the  desired  divergence  free  basic  vector  field  as 

w(x,t,cj)  =  (5^'0,-5^'0)(x,t,cj). 

A  similar  extension  can  be  constructed  for  the  scalar  field  P(x,  t,^;)  so  that  in  0o, 

P(x,t,cc;)  GL2(f^;L2(O,T;pi(0o))) 
and  becomes  —xf{t^uj)  smoothly  in  O^. 

Let  us  now  note  some  of  the  properties  of  the  constructed  basic  vector  field  w(x, 

(i)  The  flux  condition,  =  Jp  w(x,  t,^;)  •  ndP,  where  T  is  any  cross  section  of  the  channel. 

(ii)  On  the  boundary,  x  G  50,  w(x,t,C(;)  =  0,  for  all  t  G  [0,T]  and  cj  G 
(hi)  In  the  bounded  region  Oq,  the  basic  vector  field  w  satisfies 

w  G  L2(0,  T;  H^Go))  H  C([0,  T];  ^^(©o))) 

with  wt  G  L^(0,  T;  P“^(0o))).  This  gives,  in  0o 

Wf  —  uAw  +  w  •  Vw  G  L^(0,  T;  P“^(0o))).  (3.35) 

By  the  regularity  of  w  in  Oq,  we  have 

sup  |w(x, -,6^)1  =/^io(',^),  ||Vw(x,t,Cj)||L2(0,T;L2(©o))  =/^2o(^),  (3.36) 

xG©o 

where  y5io(',  •)  ^  L^(f^;  L^(0,  T))  and  /52o(')  ^  L^(f7). 

(iv)  In  the  infinite  strips  O^,  (actually  in  Oi\0o)  w{x^y,t,uj)  becomes  and  hence  satisfies 

properties  (i),  (ii)  above. 

Since,  G  ^"^([0,  T];  P^(0, 1))),  we  have 

sup  |w(x,t,w)|=  sup  \w["\y,t,uj)\  =  I3ii{uj),  (3.37) 

(x.t)eOiX(o,r)  (»/,t)e(o,i)x(o,r) 

with  pii{')  G  L^(r2). 

Also,  since  G  L^(f^;  L^(0,  T;  P^(0, 1))),  we  have 

sup  |Vw(x,  •,w)|  =  sup  \dyw["\y,-,u))\=  (3.38) 

xeOi  ?/e(o,i) 

with  p2i{’r)  eL\n;V{0,T)). 

In  O^,  is  independent  of  the  axial  direction  x,  we  have  the  constructed  vector  field  w(-,  •,  •) 
(for  almost  all  time  and  almost  all  uo  G  Cl)  has  infinite  energy  ||w||l2('0)  and  infinite  dissipation 
||Vw||l2(©)  in  0,  since 

pOO  pi  pOO  pi 

/  /  \w[^\y,t,  uj)f  dydx  =  oc  and  /  /  \dyW^\y,t,uj)fdydx  =  oo,  P  —  a.  s. 

Jo  Jo  Jo  Jo 

(v)  Also  note  that  in  (actually  in  Oi\0o), 

wt  —  uAw  +  w  •  Vw  +  VP  =  0, 


for  almost  all  a;  G  17. 
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3.1.  Connection  Between  ^  and  /3. 


Now  let  us  establish  certain  connection  between  the  norm  of  ^  and  /3  by  using  the  properties  of  w  in 
00  and  in  the  two  outlets  Oi  and  O2. 

Proposition  3.6.  /n  0  =  Oi  U  O2  U  0O;  the  regularity  properties  of  the  basie  veetor  field  w  implies  that 
/3io(-,-)  e  L2(n;L2(0,T));  /3ii(-), /3i2(-)  e  L2(a)  and  f32o{-)  e  L2(n);  ,  ■) ,  M ,  ■)  e  12(0;  1^(0,  T)). 

Thus  in  &,  we  have 


E 


[  \Plo{t,Lj)f  dt  + 
Jo 


<C(7r,z/,T)||^||2 


and 


E 


/32o{^)T  +  [  \/32i{t,u})\‘^dt+  [  \p22{t,a})\‘^  dt 
Jo  Jo 


<C{7r,i2,T)\\^\\ 


L2(il;Jfl(0,T))  ’ 


(3.39) 


(3.40) 


where  /3’s  are  defined  in  (3.36),  (3.37)  and  (3.38). 

Proof.  If  we  take  1  +  e/2  =  ?  in  the  inequality  (3.33),  one  can  see  that  for  />  <  1 

OO 

<  Ci(7r,  T,  z/)||/(-,  •)IIl2(0;L2(o,t)) 


n=0 


<C'i(7r,T,  u) 


E 

p=0 


^(2n  +  l) 

L2(0;L2(0,T)) 


4(^-6 


Now  let  us  take  the  term 
inequality  to  get 


Er=o 


L2(O;L2(0,T)) 


(3.41) 

and  use  (3.21),  (3.30)  and  Cauchy-Schwarz 


E 

p=0  *- 

=  E 


d 


L2(0;L2(0,T)) 
2 


E 

p=0 
/  OO 

=q(E 

\  \p=0 
/  OO  OO 

^MEE 

\p=0  m=0 

Coo  OO 

EE 

p=0  m=0 
/  OO  OO 

<E  EE/’^ 

\p=0  m=0 


L2(0,T)/ 


.E 

m=0 


L2(0,T)y 


L2(0,T) 


L2(0,T) 


jr--^(t,cc;) 


L2(0,T)^ 


L2(0,T) 


L2(0,T)^ 
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Coo  oo 

Eo'Eo" 

p=0  m=0 


=  E 


L2(0,T)^ 


1-p 


dt 


L2(0,T)^ 


(  '  V 

u- J 

dt 

L2(fi;L2(0,T)) 

Hence  the  inequality  (3.41)  becomes 

rT 


E 


(f 

<  Ci(7r,T,j/) 


dt 


^(2n  +  ir-6_ 

L2(0;L2(0,T))  ^  =  0 


The  series  ^  convergent  for  0  <  <^  <  |.  Hence  for  0  <  ^  <  |  and  for  p  <  1,  we  get 


E 


||^/‘fu;i(y,t,w)||^2(o,i)dtj  <  C2(7r,  T,  0  E(2n  +  1)^“^  ® 


dt 


L2(O;L2(0,T)) 


(0,T)) 

By  using  the  estimates  (3.37)  and  (3.38),  by  taking  ^  =  0,  i  in  (3.42),  one  can  easily  see  that 


(3.42) 


E 


[  \p2l{t,Lu)fdt-\-  [  \p22{t,Uj)\‘^  dt 
Jo  Jo 


=  E 


r  sup 

Jo  yeiO,l) 


\dvw\^\y,t,oj)\^  dt  + 


sup 

ye(o.i) 


IdyW^p {y,t,u})\^  dt 


<cmt)\\ 


2 

L2(r2;ffl  (0,T)) 


and 


E[{pUw)+p^M)]T 


:E 


f 


sup 

(7/,t)e(o,i)x(o,T) 


wP{y,t,u!)f  dt+  [ 
Jo 


sup  \w^^\y,t,L0)\‘^  dt 
(»/,t)e(o.i)x(o,r) 


<q|^|| 


2 

L2(r2;i^l  (0,T)) 


in  O^. 

Since  ©o  is  bounded,  by  using  property  (iii)  of  the  basic  vector  field  one  can  extend  the  constructed 
flux  ^{t^uo)  =  /p  w  •  ndS'  carried  in  into  the  domain  ©o  in  a  smooth  manner.  It  can  be  easily  shown 
that 


E  [fil^iLo)]  T  <  Cll^ll 


2 

L2(r2;ifl  (0,T)) 


and 


E 


\(3io{t,uj)\^  dt 


<C||^|| 


2 

L2(r2;ifl  (0,T))  ‘ 


Hence,  we  have  the  relation  (3.39)  and  (3.40)  in  ©.  □ 

Remark  3.7.  From  (3.38),  we  have 

sup  |Vw(x,  ■,tu)|  =  /?2i(-,u;)  G  L^(^^;L^(0,T)). 

xGOi 
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For  i  =  1,2,  we  can  always  get  an  upper  bound  for  l32i{t^uo)  for  all  t  G  [0,T]  and  almost  all  cj  G  By 
using  (3.9)  and  Young’s  inequality  for  convolution  of  two  functions,  we  obtain 


sup  sup  |Vw(x,  t,  C(;)p 

0<t<T  XGO2 


=  sup  sup 
o<t<T  ye(o,i) 


dy 


/(s,w) 


‘In  1 


00 

<  sup  47rW  / 

o<t<T 


sin(2n  +  l)7r^ 


2 


00  .  .t 

<  47r  sup  (  / 
0<t<T^  \J0 


^-2u{2n-\-lf7r^t 


2 

<  —  sup 

0<t<T 


00 


E 

n=0 


1 

(2n  +  l)2 


ds 


—  —V22{^)^ 


(3.43) 


where  rj22{^)  ^  L^(f2),  since  /  G  L^(fl;  C([0,  T]))  C  L^(fl;  L^(0,  T)). 
Hence,  we  get  sup^^[o,T]  —  V2i{^)  ^  =  1:2,  P-a.  s. 


4.  The  Linear  and  Multilinear  Operators 

In  this  section  we  define  the  Stokes  operator,  the  inertia  operator  and  the  other  operators  relevant  to  our 
analysis.  Most  of  the  results  obtained  in  this  section  has  been  taken  from  [62]  given  here  for  completeness. 


4.1.  The  Stokes  Operator 

Let  us  first  define  the  Stokes  operator  for  the  two-dimensional  admissible  channel  domain  and  analyze 
its  properties.  Let  us  denote  by  a(-,  •)  the  symmetric  bilinear  form 

a(u,  v)=  /  Vu-Vvdx.  (4.1) 

Jq 

Let  us  now  define  the  Stokes  operator  A  and  its  domain  D{A)  in  the  following  way.  Given  u  G  V,  if  there 
exists  an  element  g  G  H  such  that 

«(u,v)  =  (g,v)H,  V  V  G  V, 

then  we  say  u  G  D{A)  and  Au  =  g.  By  Poincare  inequality,  ||u||l2(©)  <  G||Vu||l2(©)  holds  for  the 
admissible  channel  domains,  we  obtain  the  coerciveness  (V-elliptic)  property  as, 

a(u,  u)  =  ||Vu||b(e)  >  C'||u||^  V  u  G  V. 

Since  the  form  a(-,  •)  is  symmetric,  continuous,  and  positive  definite,  we  have  the  following  standard 
results  obtained  in  [66]. 

Proposition  4.1.  There  exists  a  self-adjoint,  regularly  aeeretive  onto  map  A  G  Af{D{A)]'H.)  sueh  that 

a(u,  v)  =  (Hu,  v)h,  Vug  D{A),  V  v  g  V, 

with  D{A)  C  V  dense.  Moreover,  it  is  possible  to  extend  this  operator  as  an  isomorphie  onto  map 
A  G  ^(V,V')  sueh  that 

a(u,  v)  =  (Hu,  v)vx V' ,  V  u,  V  G  V. 
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The  operators  A  and  A  are  different  Friedrich’s  extensions  of  the  classical  Stokes  operator.  In  the 
remainder  of  the  paper,  both  of  these  operators  will  be  denoted  by  A.  It  follows  from  a  general  theorem 
of  Lions  [42]  that  V  =  D(A^/^).  Let  the  orthogonal  Helmhotz-Hodge  projection  Ph  :  L^(0,]R^)  ^  H 
and  Au  =  — PhAu  V  u  G  D{A). 

Let  us  now  state  a  regularity  theorem.  This  give  us  an  explicit  representation  of  D{A). 

Lemma  4.2.  Let  0  be  an  admissible  ehannel  domain.  Then,  for  a  given  g  G  R{A)  =  H,  the  Stokes  problem 
of  finding  (u,p)  :  0  ^  x  M  sueh  that 

— Au  +  Vp  =  g  m  0 
V  •  u  =  0  in  0 

=0,  u  — >  0  as  |x|  ^  oo  in  Oi  and 
j  u  •  ndP  =  0, 

has  a  unique  solution  u  that  satisfies  u  G  D{A)  =  P^(0)  H  V.  Moreover,  if  we  define  the  graph  norm 

||u||d(a)  =  ||u||l2(©)  +  ||Au||l2(©),  Vug  D{A), 
then  there  exists  C  >  0  sueh  that,  C||u||i:/2('0)  <  ||u||2:)(^)  <  ||u||j:/2('0y 

Proof  For  proof  see  Theorem  3,  section  4  of  [62].  □ 


4.2.  The  Bilinear  Operator  (Inertia  Term) 


Let  us  now  define  the  trilinear  form  and  the  bilinear  operator  and  their  properties.  The  trilinear  form  is 
given  by, 


6(u,v, 


(4.2) 


It  is  well  known  that  trilinear  form  6(',  •,  •)  :  V  x  V  x  V  ^  M  is  continuous  and  for  all  u,  v,  0  G  V  we 
have 


6(u,  V,  0)  =  — 6(u,  0,  v)  and  6(u,  v,v)  =  0. 

The  following  lemma  is  standard. 


Lemma  4.3.  [62]  There  exist  eontinuous  bilinear  operators,  P(',')  :VxV  and  P*(-,  •)  :  Vx  V  ^  V' 

sueh  that 


b{u,v,(t))  =  (B(u,v),^)v'xv  =  (•B*(v,^),u)v'xv, 


V  u,  v,(/)  G  V. 


Moreover, 

||B(u,v)||v'  <  (:^||u||L{5e)ll^"llf2je)l|v||L25e)||Vv||L^^^^^ 

||B*(v,0)||v'  <  C||Vv||L.(e)||</.||^fe)||V</.||^/4) 

and 


B*(v,</.)  =  -B*(0,v),  V</.,vGV. 


(4.3) 

(4.4) 

(4.5) 

(4.6) 


Now  we  will  provide  a  general  setting  for  estimating  the  terms  6(w,  •,  •)  and  6(-,  w,  •),  where  w  is  the 
basic  vector  field  constructed  in  Sect.  3. 


Lemma  4.4.  The  trilinear  form  6(w,  •,  •)  :  V  x  V  M  eontinuous  and  satisfies, 

(i)  6(w,  u,  v)  =  — 6(w,  V,  u),  V  u,  V  G  V  and 

(ii)  6(w,  u,  u)  =  0,  V  u  G  V. 


(4.7) 

(4.8) 
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Proof.  For  proof  see  Lemma  4,  section  4  of  [62].  □ 

Lemma  4.5.  The  form  6(-,  w,  •)  :  V  x  V  ^  M  is  continuous. 

Proof.  For  proof  see  Lemma  5,  section  4  of  [62].  □ 


Lemma  4.6.  There  exists  a  continuous  linear  operator  5i(w)  G  ^(V;  V')  such  that 

(i)  6(w,u,  v)  =  (5i(w)u,  v)v'xv,  V  u,  V  G  V. 


Moreover, 


(ii)  (Bi(w)u,v)v'xv  = -(-Bi(w)v,u)v'xv,  Vu,vgV. 
Proof.  For  proof  see  Lemma  6,  section  4  of  [62]. 


□ 


Lemma  4.7.  There  exists  a  continuous  linear  operators  52(w)  G  and  52(w)*  G  V)  such 

that 

(i)  (52(w)u,  v)v'xv  =  (u,52(w)*v)v'xv  =  6(u,  w,  v),  V  U,  V  G  V. 

Moreover, 

2 

(ii)  ||52(w)||^(V;V')  =  ||^2(w)*||^(v';V)  <  /^2o(^)  + 

i=l 

Proof.  For  proof  see  Lemma  7,  section  4  of  [62].  □ 


By  taking  the  Helmhotz-Hodge  projection  Ph,  we  get 

P(u,  v)  =  Ph(u  •  Vv),Pi(w)v  =  Ph(w  •  Vv),P2(w)v  =  Ph(v  •  Vw),  V  u,  v  G  V. 

Remark  4.8.  For  the  notational  convenience  from  now  onwards  we  use  the  symbols  fdiQ{t)  =  E[/^io(t,  cj)], 
^ij  =  E[/3ii(w)]  for  j  =  1,2,  pIo  =  E[/?|o(w)]  and  p2j(t)  =  E[P2j{t,uj)]  for  j  =  1,2. 


5.  Construction  of  the  Perturbed  Vector  Field 

Let  us  now  consider  the  system  (2.2)-(2.6)  and  introduce  the  following  change  of  variables: 

u(x,  =  v(x,t,a;)  +  w(x,  and  p(x,  t,  a;)  =  g'(x,  t,  a;)  +  P(x,  t,  a;), 

where  w(',  •,  •)  is  the  basic  vector  field.  The  regularity  properties  of  this  basic  vector  field  imply  the 
following  properties  on  Pij  and  P2j^  for  j  =  0, 1,  2: 

/32oi<^)  G  L^(a)  and  Pio{t,u)),  P2i{t,u)),  P22it,oj)  G  1^(0;  L^(0,  T)). 

Our  problem  is  to  find 

(v,  g)  ;  0  X  [0,  T]  X  a  ^  X  K 


such  that 
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Vt  +  V  •  Vv  +  V  •  Vw  +  w  •  Vv  =  —Vq  +  z/Av  +  fw  +  ^(x,  t)  in  0  x  (0,  T)  x 

V  •  V  =  0  in  0  X  (0,  T)  X 

^Ia0x[o,r]xn  =  0,  v  0  as  |x|  ^  oo  in  0^1  =  1,2, 
v(x,  0,  cj)  =  0,  (x,cj)  G  O  X  Q  and 

J  V  •  nd5'  =  0, 


where 


(5.1) 

(5.2) 

(5.3) 

(5.4) 

(5.5) 


fw  =  z^Aw  —  w  •  Vw  —  wt  —  VP.  (5.6) 

From  the  construction  of  the  basic  vector  field  it  is  clear  that  the  suppjfw}  0.  Also  one  can  notice 
that  fw  =  0  for  O^\0o.  Moreover,  form  the  estimate 

Wf  —  z/Aw  +  w  •  Vw  G  (fl;  L^(0,  T;  P“^(0o)))  , 

we  have  fw  G  (fl;  L^(0,  T;  V'))  .  A-priori  estimates  for  the  solution  can  be  obtained  by  assuming 
smoothness  of  v  and  sufficient  decay  at  infinity. 

Let  us  now  take  the  Helmhotz-Hodge  projection  on  the  Eq.  (2.7)  and  using  the  fact  that  (Vg^,  v)  =  0, 
V  V  G  H,  one  can  reduce  the  Eq.  (2.7)  to 


dv  =  (P(v(t))  +  PHfw(t))  dt  +  g{t)  dW{t), 
where  P(v(t))  =  —yAv  —  P(v,  v)  —  Pi(w)v  —  P2(w)v 
Av  =  — PhAv 
P(u,  v)  =  Ph(u  •  Vv) 

Pi(w)v  =  Ph(w  •  Vv) 

P2(w)v  =  Ph(v  •  Vw).  (5.7) 


See  Appendix  of  Mikulevicius  and  Rozovskii  [49]  for  more  details  about  the  properties  of  the  projection 
operator. 

Let  the  noise  process  be  represented  as  a  series  =  Z]/c^/c(x,t)  dll4(t,a;),  where  g  =  {gi^g2^ . . .)  is 
and  -valued  function  and  Wk  are  mutually  independent  standard  one  dimensional  Brownian  motions. 
The  stochastic  term  gdW  is  thus  an  H- valued  Wiener  process  with  a  trace-class  covariance  operator 
denoted  by  g*g  =  g*g{t)  given  by 


{9*9{t)u,v)  =  '^{9k{t),u){gk{t),v), 
k 

T^{9*9{t))  =  E  <  00. 

k 


This  means  that  the  mapping 


V 


f 


{g{t)dW{t),v) 


Vi 


{gk{t),v)  dWk{t) 


is  a  continuous  linear  functional  on  H  with  probability  1  and  the  noise  is  the  formal  time-derivative  of 
the  process  Ph^(^)  =  Jo  dW{t).  A  multiplicative  noise  of  the  form  ^(u(x,  t))  dkF(t),  where  ^(u(x,t)) 
is  a  continuous  operator  from  V  into  L^(0,  T; -^^(H)),  can  also  be  considered  as  the  random  forcing. 
Here  ^(-j  is  formally  written  and  it  is  the  time-derivative  of  g(s)  dW(s).  We  have  g(s)  dW(s)  G 
L^(f2;  C([0,  T];  H))  implies  its  time- derivative  ^(•)  satisfies  ^(•)  G  L^(f2;  W“^’^(0,  T;  H)),  since  dt  is 
linear  continuous  from  C([0,T];H)  into  T;  H). 
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6.  A-Priori  Estimates  for  the  Perturbed  Vector  Field 


Let  =  spanjei,  62, . . . ,  e^},  where  {ej}  is  any  fixed  orthonormal  basis  in  H  with  ej  G  D{A).  Let 
denote  the  orthogonal  projection  of  H  into  Define  Let  us  define  as  the  solution  of  the 

following  stochastic  differential  equation  in  the  variational  form  such  that  for  each  v  G 

d{v^{t),v)  =  {F{v^{t)),v)  dt  +  (Putwit),  v)  dt  +  '^{gk{t),v)  dWk{t), 

k 

v"(0)  =  P„v(0)  =  0,  (6.1) 

where  F{v)  =  —uAw  —  5(v,  v)  —  5i(w)v  —  52(w)v. 

It  can  be  shown  that  for  all  n  >  1,  there  exists  an  adapted  process  G  (7([0,T];  H^)  a.  s.  such  that 
satisfies  (6.1).  For  proof  see  corollary  2.1  of  Albeverio  et  al.  [1],  page  128  of  Brzezniak  et  al.  [13], 
Lemma  3.1  of  Capihski  and  Peszat  [15],  Proposition  3.2  of  Manna  et  al.  [45],  Theorem  3.1.1  of  Prevbt 
and  Rockner  [54]. 


Theorem  6.1.  Let  0  he  an  admissible  ehannel  domain.  Under  the  above  mathematieal  setting,  let  v^{t) 
be  an  adapted  proeess  in  C([0,T];  H^)  whieh  solves  the  stoehastie  ODE  (6.1).  Then  we  have  the  following 
a-priori  estimates:  For  (5  >  0 


,s)fe-^^ds 


+  ||v”(v 

<  Cl  (t,  dt^ 


and 


E 


sup  |v”(-,i)|  e 

0<t<T 


2  ^—6t 


T  2z/E 


< 


/  l|v"(' 

yo 

C2  (t,  ,  <5,  £  Tr{g*g){t)e-^*  dt 


(6.2) 


(6.3) 


Proof.  Let  us  consider  the  function  e  and  apply  Ito’s  lemma  (see  Theorem  1,  page  155  Gyongy 

and  Krylov  [26],  Metivier  [48])  to  the  process  v’^(t),  for  any  J  >  0 

d  =  -(5e-^>^pdt-2[5(v^,v^,v^)  +  6(v^,w,v^) 

+5(w,  v^,  v^)  +  z/(Av^,  v^)]  6-^'  dt  +  [(Pnfw,  v^)]  dt 

+  2  ^{gk{t),^^'^)e~^^  dWk{t)  +  e“^^Tr(^*5f)(t)  dt.  (6.4) 

k 

By  applying  the  properties  of  the  trilinear  form,  (Av’^,v’^)  =  —  and  2|(PHfw,  v’^)!  <  |^|fwp 

in  (6.4),  we  get 

d  r[e-'^‘|v”n  +  2v\W^^e-^*dt  <  dt+\\i^\^e-^^  dt 

0 

+  e~^*Tv{g* g){t)  dt  +  2  '^{gk{t),  v")e“'^*  dWk{t).  (6.5) 

k 

Let  Q  be  the  a-algebra  generated  by  in  the  probability  space  (fl;  P,  Pt,  P)-  Since  we  have  assumed  that 
the  external  random  forcing  and  the  random  flux  are  mutually  independent  processes,  the  constructed 
vector  field  w  is  independent  of  Q. 

Let  us  define  the  stopping  time  tn  by 

tn  =  inf  |f  :  |v"(f)|2  +  ||v"(s)||2  ds  >  IV 
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Note  that  tn  is  adapted  to  Q  and  w  is  independent  of  Q. 

Now  let  us  integrate  the  inequality  (6.5)  in  t  from  0  to  t  A  tat  to  get 

ptATN 

|v'"ltATA/)re  ''■'  +  -ZV  I  ds 


pi/\tn 

itATN)\^e-^^  +  2,^  Ih 

Jo 

ptATN  1 

<2j^  |(a2(wK,v")(s)|e-^Ms  +  -y^  |fwA-^*ds 

pt  ptATN 

+  /  Tr{g*g){s)e-^^ds  +  2y]  /  v-)e-^Mtyfe{s). 

Jo  .  Jo 


(6.6) 


Also  for  the  term  (52(w)v’^,  v’^)  =  6(v’^,  w,  v’^),  we  have  the  estimate 

|6(u,w,v)|  <  /32oM|u|^/2||u||^/2|v|i/^||vf/2  ^  (/32i(i,w)  +/322(t,w))  |u| 
By  using  the  above  estimate  and  using  Young’s  inequality,  we  have  for  all  G  V 

+  +  /^22(^5  ^)) 

.n\\2  ,  I  /^2o(^) 


|5(v^,av,v^)|  <  /32o(c^)|v^|||v^||  +  {(32i{t,uj)  +  (322{t,uj))  |v" 


<  V  V 


2u 


+  (J2I  +  /^22  (^5  ^) 


_n  1 2 


v"(s)||2e-‘5*  ds 


Hence  (6.6)  becomes 

ptATN 

|v”((Ar„)pe-«  +  i-  /  II' 

<  2  f 

pt  ptATN 

+  /  Tr{g*g){s)e-^^ds  +  2y'  /  (5^(5), v”)e-^MVFfe(s). 

^0  ,  Jo 


IfwPe--^*  ds 


(6.7) 


Let  us  take  conditional  expectation  on  both  sides  of  (6.7)  with  respect  to  the  cr-algebra  Q  to  get 

/  ptATN  ^ 

E  {\v^{t  ArN)fe-^^\g)  +  uE  U  ||v^(s)|pe-^^  ds  G 

P20{^) 


<  2E 


2u 


+ /32i(s,  a;)  + /322(<s,  a;)  )  ds 


+  2E 


Ifwl^e  ■’^ds 


ptATN 

Y/  (fffe(s),v")e-^«dtyfc(s) 

I,  JO 


G]+  1^  Tr{g*g){s)e-^^ds 


G\- 


(6.8) 


In  (6.8),  let  us  use  the  ^-measurability  of  and  independence  of  w  with  respect  to  G  to  get 

ptATN 


r 

\t  ATN)\^e~^^  d-u 

Jo 


n  (  ||2^— (5s 


V  s  e 


ds 


<  2E 


tATN 


(PkM 

V  2u 


+  /32i(s,w) +/322(s,w)  )  '^'*ds 


+  ^E 

0 


fwl^e  '^®ds)+  /  Tr{g*g){s)e  **  ds 


0 

ptATN 


-Ss  . 


ru/\  I  iv 

+  2yl  (5fc(s),v”)e-^*dl^fe(s). 


I.  Jo 


(6.9) 
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Let  US  take  expectation  on  both  sides  of  the  inequality  (6.9)  to  obtain 

/  ptATN 

E(|v"(i  ArAr)|^e-'5*) (  /  ||v"(s)f  e”'’*  ds 


<  2E 


+  /32i(s,  1^)  +  /322(s,  w) 


HL  V  2. 

+  Je  Ifwpe-^^  ds^  +  IV(5*fl)(s) 

/  rtATN  \ 

+  2E  f  ^  (5fc(s),v”)e-^MW^fc{s)j  . 


|v"|2e-<5« 


ds 


0} 


-'5*  ds 


(6.10) 


Note  that  the  last  term  in  the  inequality  (6.10)  is  a  martingale  with  zero  expectation,  one  obtains 

rtATN 


E 


/  f 

[|v”(tArjv)pe-^‘]  +;^E  U 


HI. 


2v 


v”(s)||2e-'^*  ds 

+  (32i{s,oj)  +  (322{s,oj)]  |v”pe“'^*ds 


fwPe-'^Ms  )  + 


fTr{g*g){s) 

Jo 


0} 


e“'^*  ds. 


(6.11) 


Now  let  us  consider 


^  rtATN  ~ 

{e( 

Jo 

/^2o(^) 

2p 


+  /32i{s,Lo)  +  /322{s,U}) 


Iv^Pe--^*  ds 


0} 


and  use  the  properties  of  conditional  expectation,  independence  of  w  with  the  a-algebra  0,  generated  by 
and  Holder’s  inequality  to  get 


rtATN  ~ 

L 

=  E  < 

l‘U’ 

=  e| 

1 

=  E<^ 

1 

f  rT 

=  E 

/  (X[ 

,Jo 

/^IqM 

2u 


+  f^21  (<5,  Uj)  +  (^22  (<5,  Uj) 
Pioi 


d5 


0} 


(X[o,tATN]|v”'(s)pe  +  (32i{s,oj)  +  (322{s,oj)  ds 


U2o(^) 

2i2 

/32o(^) 

2v 


+  p21  (<^5  ^)  +  p22  {s,  ^) 


+  p21  (<5,  Co)  +  (J22  (5,  Co) 


'1 

G  )  ds 


ds 


■ _ 2 

^  +  f^2l{^)  +  f^22{^) 


ds 


<  sup 

0<t<T 


=  sup 

0<t<T 


_ 2 

^  +  I^2li0  +  f322it) 
If  +/^2l(t)+/^22(t) 


E 


E 


(X[0.tAr«]|v"(s)|2e  dsj 

/  rtATN  \ 

^  (|v"(s)|^e-^«)  dsj. 


(6.12) 
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Note  here  that  supo<^<7^ cj)  =  rj2j{uj),  for  j  =  1,2  is  bounded  from  (3.43)  and  r]2j{’)  ^  L^(0). 
Hence,  from  inequality  (6.11),  we  get 


atATN 

||v’"(s)||2e“'5^  ds 


1, 


|fwpe-^*dsj  +  Tr{g*g){s)e-^^ds 


sup 

0<s<T 


^3^  —  — 

+  2/321(5)  +  2/322(5) 


atATN  \ 

(|v”(s)|2e-^*)  dsj  .  (6.13) 


In  particular,  we  have 


E[|v”(iATjv)|2e-^*]  <  11(3*5) (s)e-^*  ds 

atATN  \ 

(|v"(s)|2e-^«)  dsj  . 


+  sup 

0<s<T 


- _ 2 

+  2/321  (s)  +  2/^22  (^) 


An  application  of  Gronwall’s  lemma  yields 


E[K(tAriv)pe-'5*]  <  Qe  ds^  +  j\{g*g){s)e-^^  ds 


sup 

3^  —  — 

+  2/321  (t)  +  2/322  (^) 

\0<t<T 

V 

(6.14) 


By  applying  (6.14)  in  (6.13),  we  get 


X  sup 
\0<t<T 


- _ 2 

+  2/321  (^)  +  2/322  (^) 


(6.15) 


Now  let  us  take  A'  — >  oo  so  that  t  Atn  ^  t  and  by  using  the  estimates  (3.39)  and  (3.40)  to  obtain 


E  [|v”(t)|2e“^*]  +  lyE  ||v’*(s)||2e“^^  ds 

<  Cl  ||^IL.^,^,.^„,,„,5,^*Tr(3*3)(s)e-^Msj  . 


Vol.  17  (2015) 


Stochastic  Navier-Stokes  Equations  in  Unbounded  Channel  Domains 


73 


In  the  estimate  (6.9),  taking  supremum  up  to  T  A  tn  before  taking  the  expectation,  we  obtain 


E 


sup  |v"(t)|2e-^* 

jd<t<T  Atn 


+  j/E 


,n (+\  \\2^  —  6t 


e-'’*  dt 


j-  l  ■/\TN 

/  " 

Jo 

[  ^(flfe(s),v”(s))e“'^^dWfe(s) 

Jo  , 


dt 


sup 
0<t<T Atn  • 


t  /  o2 


/^2o(^) 


I  sup  ^(  [  (  ^ 

[0<t<TATN  \Jo  \ 


-/?2l(s,tu)  +/?22(s,t^)  1  |v"|^e 


(6.16) 


Now  by  applying  Burkholder-Davis-Gundy  inequality  and  Young’s  inequality  to  the  term 
2E  [supo<t<rATjv  lo  Efe(flfc(s)>  v”(s))e-'5^  dlCfe(s)]  ,  one  gets 


2E 


sup 

0<t<TATN 


< 


[  Y.^9k{s),v-{s))e-^^dWk{s) 

■^0  k 

(pT At n  \ 

Tr(5*5)(t)|v"(t)|2e-2^Mt 


1/2 


<  2V2E 


/  \  \ 
sup  |v”(t)|2e-^M  /  Tr{g*g){t)e-^^dt] 

\0<t<TATN  /  \Jo  J 


1/2- 


<iEf  sup  \v^{t)\^e-^*]  +4  [  Tr{g*g){t)e-^*  dt. 
2  \0<t<TArN  )  Jo 


(6.17) 


From  (6.12),  it  can  be  easily  seen  that 


E I  sup  ^  (  [ 

lo<t<TArAr  \Jo 


t  r  o2 


/^2q(^) 

2ia 


+  /^21  (<^5  ^)  +  /^22  (<5,  ^) 


d5 


<  sup 

0<t<T 


- _ 2 

^  +  /^22(^) 


(r 


•TAtat  \ 

E|  /  sup  ^"(5)126-5^5  . 

0<s<t  / 


Hence  (6.16)  becomes 


(6.18) 


E 


sup  |v"(t)|2e-^* 

_0<t<TATN 


2z/E 


<  -E 

^  V  ./o 


f 


IfwPe  1  +  10 


2  sup 

0<t<T 


^  —  — 

^+fJ2l{t)+P22{t) 


[L 

rTr{g*g){t) 

Jo 

I 


dt 


dt 


TAtn  ^ 

E|  /  sup  |v"(s)|A-'^®ds 

0  0<s<t  j 


(6.19) 
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In  particular,  we  have 
E 


sup  |v”(t)|2e-'^* 
_0<t<T  Atn 


2^ 


{I. 


T  \  fT 

2„-St 


|fwre-^*di  +10 


/  Tv{g*g){t)i 
Jo 


e-^*  dt 


sup 

0<t<T 


3^  —  — 

+  2/^21  (t)  +  2/^22  (^) 


E 


/  pTAtn  / 

\Jo  V' 


sup  |v^(5)|^e-'^"  )  d5 

0<s<t 


Applying  Gronwall’s  inequality,  one  gets 


E 


sup  |v"(i)|2e-^* 

P<t<TArN 


{i: 


<  I  9E  I  /  +10 


rTr{g*g)it) 

Jo 


dt 


X  exp  sup 

\0<t<T 


_ 2 

+  2/321  (t)  +  2/322  (^) 


Hence  by  applying  (6.20)  in  (6.19),  we  deduce  that 


E 


pTAtn 


nl  /\Tn 

sup  |v^(t)pe-^^  +2z/  /  E||v^(t)|pe-^Mt 

_0<t<T  Atn  _  -'O 


e-'5*  dt 


1  +  2  exp  sup 

\0<t<T 


—  — 

+  2/321  (t)  +  2/322  (^) 


X  I  sup 
\  0<t<T 


6^  —  — 

+  2/321  (t)  +  2/322  (^) 


Let  us  now  take  N  ^  oo  so  that  T  A  tn  ^  T  and  by  using  (3.39)  and  (3.40)  to  get 

rT 


E 


sup  Ivl^e-^' 

0<t<T 


■  2z/E 


/  l|v||^ 
yo 

<C2  . 


(6.20) 


(6.21) 


Hence  we  get  the  required  result. 


□ 


Theorem  6.2.  Let  0  be  an  admissible  ehannel  domain.  Under  the  above  mathematieal  setting,  let  w^{t) 
be  an  adapted  proeess  in  C([0,T];  H^)  whieh  solves  the  stoehastie  ODE  (6.1).  Then  we  have  the  following 
a-priori  estimates: 


E|v"(.,t)+i/E  y*||v-(.,s)||2ds) 

<C3  (6.22) 


and 
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e(  sup  + 

\0<t<T  Jo  j 

<  C4  .  (6.23) 

Proof.  Apply  Ito’s  lemma  (see  Theorem  1,  page  155  Gyongy  and  Krylov  [26],  Metivier  [48])  to  the  function 
and  to  the  process  v’^(t).  Proceeding  similarly  as  in  Theorem  6.1  and  using  the  inequality 

2|(PHfw,v")|<|K||2  +  ||f^||^„ 

we  get  the  required  result.  □ 


7.  Existence  and  Uniqueness  of  Strong  Solutions  for  the  Perturbed  Vector  Field 

Monotonicity  arguments  were  first  used  by  Krylov  and  Rozovskii  [34]  to  prove  the  existence  and  unique¬ 
ness  of  the  strong  solutions  for  a  wide  class  of  stochastic  evolution  equations  (under  certain  assumptions 
on  the  drift  and  diffusion  coefficients),  which  in  fact  is  the  refinement  of  the  previous  results  by  Pardoux 
[50,51]  (also  see  Metivier  [48])  and  also  the  generalization  of  the  results  by  Bensoussan  and  Temam  [10]. 
Menaldi  and  Sritharan  [47]  further  developed  this  theory  for  the  case  when  the  sum  of  the  linear  and 
nonlinear  operators  are  locally  monotone. 

In  this  section,  we  will  prove  the  local  monotonicity  of  the  sum  of  the  Stokes  operator  and  the  inertia 
term  of  the  perturbed  vector  field,  and  following  Menaldi  and  Sritharan  [47]  we  use  a  generalization  of 
Minty-Browder  technique  to  prove  the  existence  and  uniqueness  result  that  avoids  compactness  method 
and  hence  applicable  for  unbounded  domains  directly. 

We  start  this  section  with  a  lemma  known  as  the  Gagliardo-Nirenberg  inequality  (section  1.2,  Theorem 
2.1  in  DiBenedetto  [20]). 

Lemma  7.1.  Let  cp  G  1T^’^(M’^),  where  n  is  the  dimension  of  the  spaee.  For  every  fixed  number  p,  s  >  1, 
there  exists  a  eonstant  C  depending  only  upon  n^p  and  s  sueh  that 

where  a  G  [0, 1],  p,  g'  >  1  and  s  satisfies  the  following  relation: 


The  following  lemma  is  a  special  case  of  Lemma  7.1  in  two  dimension  which  will  be  useful  in  our 
context. 

Lemma  7.2.  For  p  G  (7^(0),  where  0  is  the  admissible  ehannel  domain,  we  have  the  following  estimate: 

||p||l4(©)  <  2||(p||l2(©) II V(p||l2(©). 

The  general  proof  in  of  this  lemma  can  be  obtained  from  Ladyzhenskaya  [36]  (Ghapter  1,  Lemma  1). 

Remark  7.3.  The  above  lemma  suggests  that  V  D  H  C  L4(e)  and 

L2(0,  T;  V)  n  L~(0,  T;  H)  C  L^((0,  T)  x  0).  (7.1) 

Definition  7.4.  Let  X  be  a  Banach  space  and  let  X  be  its  topological  dual.  An  operator  F  :  D  ^  X  ,D  C 
X  is  said  to  be  monotone  if 

{F{x)  -  F{y),  X  -  y)^/  xX  ^  0  x,y  e  D. 

F  is  said  to  be  A-monotone  if  F  +  A/  is  monotone,  where  A  G  M  and  /  is  the  identity  operator. 

Next  we  prove  the  local  monotonicity  property. 
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Lemma  7.5.  (Local  Monotonicity  of  F)  Let  us  denote  by  the  elosed  L^-ball  in  V; 

=  {x  G  V  :  ||x||l4  <  g}. 

Define  the  nonlinear  operator  F  on  V  by  F{v)  =  —uAv  —  —  B2{'w)v  —  .  Then  the  operator 

F  +  KI,  for  A  G  M  and  suffieiently  large  |A|^  is  monotone  in  Mg,  i.e.,  for  any  v  G  V  and  x  G  ®^, 

(F(v)-F(x),z)  +  A|z|^  <0,  (7.2) 

where  z  =  v  —  x  and  A  =  ^  —  —  |r?2i(w)|  ~  |??22(<^)|  “  <^20  is  defined  in  (3.36)  and  ri2i ’s  are 

defined  in  Remark  3.7. 


Proof.  We  have  F{v)  =  —uAv  —  B(y,  v)  —  52(w)v  —  5i(w)v  and  F(x)  =  — nAx  —  5(x,  x)  —  52(w)x  — 
5l(w)x. 

By  a  simple  application  of  the  property  of  the  trilinear  form,  we  have 

(52(w)v  -  52(w)x,  Z)  =  (52(w)v,z)  -  (52(w)x,  z) 

=  6(v,  w,  z)  —  6(x,  w,  z) 

=  6(v  —  X,  w,  V  —  x)  =  6(z,  w,  z) 

=  (52(w)z,z). 

Similarly, 


(5i(w)v  —  5i(w)x,  z)  =  (5i(w)z,z). 


Now  one  can  get 

(F(v)  -  F(x),z)  =  (z/Az,z)  -  (B(v,v)  -  B(x,x),z)  -  (B2(w)z,z)  -  (Bi(w)z,z). 

Also  we  have  (z/Az,z)  =  — iv||z||^  and 

(B(v,z),x)  =  6(v,z,x)  =  -6(v,x,z)  =  -6(v,x,z)  -  6(v,z,z) 

=  -6(v,  X  +  z,  z)  =  -6(v,  V,  z) 

=  -(•B(v),z). 

Similarly,  one  can  prove  that 

(a(v),z)  = -(B(v,z),x)  and  (a(x),  z)  = -(B(x,  z),  v). 


(7.3) 


This  gives 

(a(v)  -  B(x),  z)  =  -(a(v,  z),  x)  +  (B(x,  z),  v)  =  6(x,  z,  v)  -  6(v,  z,  x) 
=  6(x,  z,  v)  —  6(x,  z,  x)  +  6(x,  z,  x)  —  6(v,  z,  x) 

=  6(x,  z,  V  —  x)  —  b(y  —  X,  z,  x)  =  6(x,  z,  z)  —  6(z,  z,  x) 
=  -6(z,z,x)  =  -(B(z),x). 

Using  Holder’s  inequality  and  Sobolev  embedding  theorem,  one  gets 

|(B(v)  -  B(x),z)|  =  I  -  (B(z),x)|  =  |6(z,z,x)| 

<  l|z||L4(e)||z||v||x||L4(e) 

<  ||zf/2|zy2||x||j^4^Q^ 

Also,  we  have  (Hi(w)z,z)  =  6(w,z,z)  =  0.  By  using  the  inequality 

|6(u,W,v)|</320|ur/2||ur/2|v|V2||v||l/2  +  ^^^^(^)|u||^|^ 
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Young’s  inequality  and  Remark  3.7,  we  get 


|(-B2(w)z,z)|  =  |&(z,  w,z 


<  /?2o(t^)|z|||z||  +  '^/32k{t,0j)\ 

/32o(^) 


*=1 


+  \r]2iioj)  \  +  \r]22{oj)\ 


Applying  all  these  estimates  in  (7.3),  one  can  deduce  that 

(F(v)  -  F(x),z)  <  -|||z||2  +  +  \ri2i{uj)\  +  |??22H|^  |zp- 

Since  0  is  an  admissible  channel  domain  and  Poincare  inequality  still  holds,  we  get  |z|  <  C||z||  and  hence 
we  have 

(F(v)  -F(x),z)  +  -  \ri2i{uj)\ -  |??22HI  -  |zp  <  0.  (7.4) 

□ 


Definition  7.6.  {Strong  Solution)  A  strong  solution  v  is  defined  on  a  given  probability  space 

as  a  L^(0,  T,  V)  H  C(0,T;H))  valued  adapted  process  which  satisfies  the  stochastic  channel  flow 

model 


dv  +  [5(v,  v)  +  52(w)v  +  5i(w)v]  dt  =  [-vAv  +  Pr^w]  dt  +  g{t)  dW{t) 

v(0)  =  0, 


in  the  weak  sense  and  also  the  energy  inequalities 


e(  sup  \v{-,t)\^ +2v  [  ||v(-,i)||2ds  ) 

\0<t<T  Jo  J 


and 


E 


sup  |v(-,i)|2e-'5* 

0<t<T 


<  c 


T,y,\\-n 


L2(r2;7fl(0,T)) 


Jo 


Tr(5*5)(^)e 


-5t 


(7.5) 


(7.6) 


(7.7) 


Definition  7.7.  Let  X  be  a  Banach  space  and  let  X'  be  its  topological  dual.  An  operator  F  :  D  ^  X' ,  D  C 
X  is  said  to  be  hemicontinuous  at  x  G  P,  if  ^  G  X,  >  0,  n  =  1,  2, . . .,  — >  0  and  x  +  tny  G  D  imply 

F{x  F  tnV)  F{x)  weakly. 

Theorem  7.8.  Let  fw  G  L^(0,T;  V').  Then  there  exists  a  unique  adapted  proeess  v(x,  Lcj)  with  the  regu¬ 
larity 

V  G  C(0,  T;  H)  n  E(0,  T;  V)) 

satisfying  the  stoehastie  model  (7.5)  and  the  a-priori  bounds  in  (7.6)  and  (7.7). 


Proof.  Part  I:  Existence  of  strong  solution 

We  prove  the  existence  of  Strong  solutions  of  the  stochastic  channel  flow  model  (7.5)  in  the  following 
four  steps. 

Step  (1)  Finite- dimensional  (Galerkin)  approximation  of  the  stoehastie  ehannel  flow  model  (7.5) 


78 


U.  Manna  et  al. 


JMFM 


Let  {ei,  62, . . . ,  e^, . . .}  be  a  complete  orthonormal  system  in  H  belonging  to  V.  Denote  by  the 
n- dimensional  subspace  of  H  generated  with  {ei,  62, . . . ,  e^}.  Let  be  the  solution  of  the  following 
stochastic  differential  equation  in 

d{v^{t),v)  =  {F{-v^{t)),v)  dt  +  (PHfw(i),  v)  dt  +  '^{gk{t),v)  dWk{t), 

k 

v”(0)  =  P„v(0)  =  0,  (7.8) 


in  (0,T)  for  any  v  G  where  F(y)  =  —uAv  —  5(v,v)  —  5i(w)v  —  52(w)v.  Denoting  F{v)  = 
F(v)  +  ^Hfw,  we  have  {v’^}  satisfies  the  stochastic  Itb  differential  equation 

dw^{t)=F{w^{t))dt^g{t)dW{t),  v^(0)  =  0  (7.9) 

and  the  corresponding  energy  equality 

d|v”(i)  |2  =  2(F(v"(i)),  v”(i))  dt  +  2  {g{t)  dW{t),  v"(t))  +  Tr{g*g){t)  dt.  (7.10) 

Step  (2)  Weak  convergent  sequences 

Using  the  a-priori  estimates  in  Theorems  6.1  and  6.2,  it  follows  from  the  Banach-Alaoglu  theorem 
that  along  a  subsequence,  the  Galerkin  approximations  {v’^}  have  the  following  limits: 

v”  ^  V  weak  star  in  L2(11;  L°°(0,  T;  H)  n  L2(0,  T;  V)) 

F(v")  ^  Fo  weakly  in  L'^{n;  1^(0,  T;  V')).  (7.11) 

Now  the  assertion  in  the  second  statement  holds  since  F(v")  is  bounded  in  the  space  L^(n;  L^(0,T;  V')). 
Note  that  F(v”(t))  ^  Fo(t)  =  Fo(t)+PHfw(t)  weakly  in  (11;  L2(0,T;V'))  since  fweL2(ll;L2(0,T;V')). 

Moreover,  v”(t)  €  L^(n;  (^([OjT];  H„))  and  the  Galerkin  approximations  v”  converge  to  v  weakly  star 
in  the  Banach  space  L^(n;  C'([0,  T];  H))  implies  that  t  i— >  v(t)  is  a  continuous  function  from  [0,T]  into 
H  with  probability  1  (see  page  46-47,  Proposition  3.3  of  Menaldi  and  Sritharan  [47]).  Hence  we  have 
v(t)  G  L2(11;  C(0,  T;  H)  n  L2(0,  T;  V)). 

Now  let  us  apply  Ito’s  lemma  to  the  function  and  to  the  process  v”'(t)  to  obtain 

d[e-’-d)|v"(t)|2]  =  -e-’’(*)[(r(t)v"(t),v"(t))  +2(P(v"(t)),v"(t))]  dt 

+  2e-^^*\g{t)  dlP(t),  v”(t))  +  e- ^ Tr {g* g){t)  dt,  (7.12) 


where  r{t)  denotes  the  derivative  of  r{t).  Let  us  integrate  the  above  equality  (7.12)  from  0  to  T  and 
taking  the  expectation  to  get 


E 


(T)l^ 


=  E 


J  (2F(v”(t))  -  r(t)v”(t),  v”(t))  dt 

Te-’'W(5(t)dlP(t),v”(t))  +E  Te-''WTV(5*fl)(t)dt 

Jo  Jo 


But  2  e  idit)  blU (t),  v’^(t))  is  a  martingale  having  expectation  zero.  Then  from  the  above  equation, 
we  have 


E 


=  E 


^  (2FO 


‘(t))-r(t)v”(t),v^ 


dt 


+  f  e  '^^*'’Tr{g*g){t)dt 

Jo 


From  (7.11)  and  (7.9),  we  also  note  that  v  has  the  Ito  differential 

dv(t)  =Po(i)dt  +  5(t)dlP(t)  weakly  in  1^(11;  1^(0, T;V')), 


(7.13) 


(7.14) 
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with  v(0)  =  0,  where  Fo(^)  =  ^o{t)  +  Also  v  satisfies  the  energy  equality 

d|v(t)|^  =  2{FQ{t),v{t))  dt  +  2  {g{t)  dW{t),v{t))  +Tv{g*g){t)dt.  (7.15) 


By  applying  Ito’s  lemma  to  the  function  e  to  the  process  v{x,t,u>)  from  (7.14),  we  can  show 

that  satisfies 


E 


.-riT) 


v(T)|' 


=  E 


f 


2Fo{t)  -  r(t)v(t),  v(f) )  dt 


+  /  e-^^*^Tr{g*g){t)dt 
Jo 


(7.16) 


Step  (3)  Lower  semicontinuity  property 

Since  — >  v  weak  star  in  L^(0,  T;  H)  D  L^(0,  T;  V)),  by  the  lower  semicontinuity  property  of 

L^-norm,  we  have 


lim  inf  E 


,-r(r)|^n(j.)|2 


>  E 


;-»’(r)|v(r)|2 


Let  US  take  lim  inf  on  the  Eq.  (7.13)  and  by  using  (7.17)  and  (7.16)  to  obtain 

n 

jT  e-r(t)  (2F(v’^(f)) -r(f)v’^(f),v’^(f))  df 


Hence,  we  get 


lim  inf  E 


=  lim  inf  E 


o-riT)\ 


>  E 


=  E 


v”(T)|2-  [  e-''WTf(ff*ff)(t)df 

Jo 

(T)|2-  r  e--^*'>Tv{g*g){t)dt 
Jo 

J  i^Foit)  -  r(f)v(t),  v(f))  df 


lim  inf  E 


>  E 


jT  e-Ht)  (2F(v"(f))  -  r(f)v"(f),v' 
(2Fo{t)-r{t>{t),v{t))  dt 


dt 


Step  (4)  Local  Minty-Brow der  technique 

For  y(cc;,x,  t)  G  L^(fl;  L'^(0,  T;  H^))  with  m  <  n,  let  us  set 


r{t)  =  r{t,w)  =  2C{u})t  +  ^  ^(dllb  ds, 


(7.17) 


(7.18) 


(7.19) 


(7.20) 


where  C{lo)  =  +  \  V2i{^)  \  +  \  V22{^)  \  —  ^  as  an  adapted,  continuous  (and  bounded  in  uj)  real- valued 

process  in  [0,T].  Note  that  from  Lemma  7.2  and  Remark  7.3,  r{t^uo)  is  well  defined. 

Since  y(t)  =  y(a;,x,  t)  G  L^(fl;  L^(0,  T;  H^)  belongs  to  the  closed  L^-ball  in  V,  from  the  Local 
Monotonicity  Lemma  7.5,  by  setting  z(t)  =  v’^(t)  —  y(t)  and  F(y^)  =  F(v’^)  +  PRfw  with  fw  G 
L^(0,T;  V'),  we  have 
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E 


(2F(v"(t))  -  2F{y{t)),v^{t)-y{t))  dt 
r  e-^^*^r{t)\v^{t)  -  y{t)\^  dt 

Jo 


-E 


On  rearranging  the  terms  in  the  inequality  (7.21),  we  get 

pT 


<  0. 


E 


^  (2F(v”(t))  -  r(t)v"(t),  v”(t)  -  y(t))  dt 

(^2F{y{t))  -  r(t)y(t),  v"(i)  -  y{t))  dt 


Taking  the  limit  as  n  ^  oc  in  (7.22),  one  obtains 

pT 


lim  inf  E 


<  E 


^  (2F{v^{t))  -  f{t)v-{t),v^{t)  -  y{t))  dt 

{2F{y{t))  -  r{t)y{t)Mt)  -  y{t))  di 


By  using  (7.19)  and  (7.23),  we  get 

pT 


E 


<  E 


From  the  above  inequality,  we  get 

pT 


^  e-’-d)  (2Fo{t)  -  r{tMt),  v(t)  -  y(t))  dt 

jT  e-r-m  (^2F(v”(t))  -  r(t)v”(t),v”(t)  -  y(t))  dt 
(2F(y(t))  -  r(t)y(t),  v(t)  -  y(t))  dt 


<  lim  inf  E 


E 


^  e-'-d)  (2Fo(t)  -  r(t)v(t),  v(t)  -  y(t))  dt 

£  (2F(y(t))  -  r(t)y(t),  v(t)  -  y(t))  dt 


<  E 


On  rearranging  the  terms  in  the  inequality  (7.25),  we  obtain 

pT 


E 


^  e-'-d)  (2Fo{t)  -  2F(y(t)),v(t)  -  y(t))  dt 

re-’-Wr(t)|v(t)-y(t)|2dt 

Jo 


(7.21) 


(7.22) 


(7.23) 


(7.24) 


(7.25) 


(7.26) 


This  estimate  holds  for  any  y  G  L^(r2;  L°°(0,  T;  H^))  for  any  m.  G  N.  It  is  clear  by  a  density  argument 
that  the  above  inequality  remains  the  same  for  any  y  G  L^(f2;  C(0,  T;  H)  nL^(0,T;V)).  Indeed,  for  any 
y  G  L^(fl;  C(0,  T;  H)  H  L^(0,  T;  V)),  there  exits  a  strongly  convergent  sequence  y^  G  L^(fl;  C(0,  T;  H)  H 
L^(0,T;  V))  that  satisfies  inequality  (7.26). 

Let  us  now  take  for  any  A  >  0,  y(t)  =  v(t)  —  Az(t)  in  (7.26).  Note  that  v  G  L^(fl;  C(0,  T;  H)  D 
L^(0,T;V))  satisfies  the  ltd  differential  equation  given  by  (7.14)  and  z  is  an  adapted  process  in 
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L2(n;C(0,T;H)  nL2(0,r;V)).  Then  from  (7.26),  we  have 


AE 


2Fo{t)  —  2F(y(t)  —  Az(t)),  z(t) )  dt 


-U 


o-rit) 


mmfdt 


<  0. 


Dividing  by  A  on  both  sides  of  the  above  inequality  and  letting  A  go  to  0  and  using  the  hemicontinuity 
property  of  F,  one  obtains 


E 


[Poit)  -  Fivit)),zit)) 


dt 


<  0. 


(For  more  details  see  Fernando,  Sritharan  and  Xu  [22].)  Since  z  is  arbitrary,  we  conclude  that  Fo(t)  = 
F(v(t)).  Thus  the  existence  of  the  strong  solutions  of  the  stochastic  system  (7.5)  has  been  proved. 


Part  II:  Pathwise  uniqueness  of  strong  solution 

Let  vi  G  L^(r2;  C(0,  T;  H)  H  L^(0,  T;  V))  be  another  solution  of  the  Eq.  (7.5).  Then 'd  =  v  —  vi  solves  the 
stochastic  differential  equation  in  L^(f];  L^(0,  T,  V')),  given  by 

d^  =  (F(v(t))-F(vi(t)))dt,  ^(0)  =  0.  (7.27) 

On  multiplying  (7.27)  by  we  obtain 

d[e-"(*VWI^]  =  -e-"(*V(f)|t?(f)|2df  +  2e-’’(*)(F(v(f))  -F(vi(f)),i9(f))dh 

Now  let  us  use  the  local  monotonicity  condition  (7.2).  Remark  7.3  gives  Vi  in  the  L^-ball  of  V,  then  by 
taking  z  =  v  —  vi  in  (7.4)  and  by  choosing  r{t)  as  in  (7.20),  one  gets 

d[e-’’(*V(t)|2]  <  0  ^  e-’'(*V(t)|2  <  |i?(0)|2  =  0. 

Hence  the  uniqueness  of  v(x,t,C(;)  satisfying  the  system  of  Eqs.  (2.7)-(5.6)  has  been  proved.  □ 


8.  Characterization  of  the  Perturbation  Pressure 

In  this  section,  we  characterize  the  perturbation  pressure,  g(',  *,  •)  =  p(-,  •,  •)  —  P(',  •,  •)  by  using  a  general¬ 
ization  of  de  Rham’s  theorem  [56]  to  processes  (see  [40]  for  more  details).  We  have  constructed  the  basic 
vector  field  in  such  a  way  that  the  random  basic  pressure  field  P(x,  t,  cj)  tends  to  infinity  in  each  of  the 
outlets  as  |x|  ^  oo. 

Theorem  8.1.  [A  generalization  of  de  Rham’s  theorem  to  processes]  Let  D  be  a  hounded,  eonneeted  and 
Lipsehitz  open  subset  o/M^.  Given  ro  G  [l,oo],  ri  G  [l,oo]  and  si  G  Z,  let 

X  e  (fl; W"^’^i(0,T;P-^(P)))  (8.28) 

be  sueh  that,  for  all  v  G  ^{D),  V  •  i;  =  0,  F-a.  s., 

(X5'^)^'(d)x^(d)  =  0  m  ^'(0,T),  (8.29) 

where  ^'{D)  denotes  the  dual  spaee  of  S^{D).  Then,  there  exists  a  unique 

q  G  (fl;  W"^’^^(0,T;L2(P)))  (8.30) 

sueh  that  F-a.  s., 

=  X  m  ^'((0,T)  X  P),  (8.31) 

I  g'dx  =  0  m  ^'(0,T).  (8.32) 

Jd 

Proof.  Eor  proof  see  Theorem  4.1  of  [40].  □ 
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(8.33) 

(8.34) 


Lemma  8.2.  Let  0  be  an  admissible  channel  domain,  letQcQ  (0  is  non-empty)  be  a  bounded  subdomain 
with  0  C  0.  For  any  given  tq  G  [1,  oo],  ri  G  [1,  oo]  and  si  G  Z,  let 

XGL"°  (a;W*'’’'HO,T;F-i(0))) 
be  such  that,  for  all  v  G  ^(0),  V  •  'i;  =  0;  F-a.  s., 

(X:'^)^'(©)x^(©)  =  0  m  ^'(0,T). 

Then  there  exists  a  unique 

gGL’-o(a;W*^’’-K0,T;LL(e))), 

satisfying 

X  =  Vg 

in  the  sense  of  distributions  in  (0,  T)  x  0,  and 


L 


gdx  =  0 


for  almost  all  t  G  [0,T]. 


The  proof  for  deterministic  case  can  be  found  in  Lemma  1.4.2,  Chapter  2,  page  202  of  Sohr  [58]  for 
any  general  unbounded  domain  >  2.  Also  see  the  Remark  4.3  of  Langa  et  al.  [40]  for  stochastic 

case. 

We  will  now  prove  the  existence  and  uniqueness  of  the  perturbation  pressure  using  Lemma  8.2. 
Theorem  8.3.  There  exists  a  unique  scalar  distribution 

qeL^  (n;W-i’~(0,T;LL(e))) 

such  that  (2.7)-(5.6)  are  satisfied  in  the  sense  of  distributions. 

Proof.  By  the  existence  and  uniqueness  theorem  on  strong  solutions,  there  exists  a  process  v{-K,t,uj) 
satisfying  (2.7)-(5.6)  and 

v{^,t,Lu)  G  L^  (fl;C(0,T;H)  nL2(0,T;  V))  .  (8.35) 

Also  v{x,t,uj)  satisfies  the  variational  equation  P-a.  s.,  for  all  t  G  [0,T] 

•Zdx  =  — z/  /  f  Vv(s)  •  V^^dxds  —  [  /  (vfs)  •  Vvfs))  •  Z  dxds 

/©  Jo  Jq 

pt  p  nt 

•Zdxds—  /  /  (wfs)  •  Vvfs))  •  Z  dxds 


[  v(t)-Zdx  =  — z/  f  f  Vv(s)  •  V^^dxds  —  f  f  {v{s)  •  \/v{s)) 

Je  Jo  Jq  Jo  Jq 

-  f  f  (v(5)  •  Vw(5))  •  Zdxds  -  f  f  (w(5)  •  Vv(5)) 
Jo  Jq  Jo  Jq 

+  [  (fw(5),  ^)i^-i(©)xi^i(©)  ds  +  [  {g{s)dW{s),Z), 
Jo  Jo 


(8.36) 


for  all  Z  G  Hq{&)  such  that  V  '  Z  =  0.  Let  us  differentiate  (8.36)  with  respect  to  t  {uj  Q  ft  being  fixed), 
we  get  in  ^'(0,  T) 


—  /  dtv-Zd:si  —  u  I  Vv-V2^dx—  /  (v-Vv)-Zdx—  /  (w-Vv)-Zdx 

Jq  Jq  Jq  Jq 


-  /  {-v  -S/w)  •  Zdy^P  {^^{J,Z)H-, 


iQ)xH^iQ) 


ue)P{^{J.Z)=0. 


Since  Z  G  i7o(0),  and  (Vv,  V2^)  =  —  (Av,Z) ,  from  the  Eq.  (8.37),  we  have 

(-5tv  +  z/Av  -  V  •  Vv  -  V  •  Vav  -  Av  •  Vv  +  fw  +  ^(x,  t),  Z) =  0. 

Let  us  denote  y  =  —dfV  +  z^Av  —  v  •  Vv  —  v  •  Vw  —  w  •  Vv  +  +  ^(x,  t).  We  will  prove  that 

XGL'  (f2;W-i’-(0,T;  77-1(0))). 


(8.37) 

(8.38) 


(8.39) 
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Since  is  linear  continuous  from  L^(0, T;  H)  into  W  ^’^(0, T;  H)  and  then  into  W  ^(0)), 

by  using  (8.35),  we  have 

dtv  e  W-^’^(O,T;iJ-\0)))  . 

We  know  that  A  is  a  linear  continuous  operator  from  i^^(0),  and  then  from  V,  into  Hence 

(8.35)  implies  that 

z/Av  G  {n;L\0,T,H-\e)))  C  (f^;  L1(0,  T,  i^-\0)))  . 

Now  for  every  v  G  L^(0,  T;  i^“^(0)),  we  have  v  =  dt  fgV  and  /q  v  G  L^(0,  T;  i7“^(0)).  This  gives  v  G 
pp-i,oo(o  7";  ^-1(0))^  by  using  the  definition  of  ^“^’"^(0,  T;  i^“^(0)).  Hence,  we  have  the  topological 
embedding 

L^(O,T;i7-^(0))  C  W"^’^(O,T;i^-\0)).  (8.40) 

Hence,  we  get 

i/Av  e  (a;Wi’~(O,T;iJ-i(0)))  . 

From  Lemmas  4.3,  4.4,  4.5  and  4.7,  it  is  straight  forward  that  the  terms  v  •  Vv,  v  •  Vw  and  w  •  Vv 
are  in  (f];  L^(0,  T;  i7“^(0)))  which  is  clearly  included  in  the  space  L^(f^;  W“^’'^(0,  T;  i7“^(0))). 

The  term  =  0  in  Oi\0o  and  from  the  estimate  Wf  —  z/Aw  +  w  •  Vw  G  (fl;  L^(0,  T;  L^(0o)))  C 
(f);L2(O,T;i7-i(0o)))  ,  we  have 

fw  G  L2  (Q-,L^(0,T-,H-^(e)))  C  (fl-,L\0,T-,H-^(e))) . 

Hence,  by  the  embedding  (8.40),  we  get  fw  C  (^4;  W“^’^(0,T;  i7“^(0))). 

Since  jQ^(5)dW(s)  G  L^(r2;  C([0,  T];  L^(0))),  its  time-derivative  ^(•,  •,  •)  satisfies 
^(x,  t,a;)  G  (O;  W“^’^(0,  T;  L^(0)))  ,  which  is  clearly  included  in  the  required  space.  Hence,  we  have 
(8.39). 

Finally  by  applying  Lemma  8.2,  we  get  a  unique 

geLi  (n;W-i’~(0,T;LL(e))) 

such  that  X  —  ^Q  fhe  Navier  Stokes  equations  are  satisfied  in  the  distributional  sense.  □ 


9.  Conclusion 

Theorem  9.1.  Let  0  be  an  admissible  ehannel  domain.  Then  for  any  given  random  net  flux  ^{t^uj)  G 
L^(f2;  iL^(0,  T))  through  the  ehannel  domain,  there  exits  a  unique  pathwise  strong  solution  u(x,  t,^;)  and 
pressure  p{yi,  t,  uo)  sueh  that, 

u  -  w  G  L^  (f];  C(0,  T;  H)  H  L^{0,  T;  V)  n  L^(0,  T  :  H)) 

and 


P-PgL^  (n;W-i’~(0,T;LL(e))) 

of  the  problem  (2.2)-(2.6);  where  w  is  the  eonstrueted  basie  veetor  field  with  Jp  u  •  ndS'  =  Jp  w  •  ndS'  = 
^{t,uj)  for  eaeh  uj  E  ft  and  P  is  the  eonstrueted  sealar  field. 

Proof.  The  existence  of  the  solution  u(x,  t,  cj),  satisfying  the  prescribed  flux  condition  ^{t,  u)  is  given  by 
the  existence  of  the  basic  vector  field  w(x,t,a;),  which  carries  the  flux  ^{t,u))  and  the  divergence  free, 
zero  flux  vector  field  v(x,  t,^;). 
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For  the  pathwise  uniqueness  of  the  solutions,  let  us  assume  that  there  exists  two  solutions  ui(x,  t,^;) 
and  U2(x, satisfying  the  prescribed  flux  condition  Hence  the  difference  field  tv{t)  = 

m(x,  =  ui(x,  —  U2(x,  t,^;)  satisfies  the  equation 


dm(t)  =  —iyAtv{t)  dt  —  [H(m(t),  tu(t))  +  H(ui(t),  tu(t))  +  U2(t))]  dt, 

m(0)  =  0, 

to  •  nd5'  =  0,  V  cj  G 


(9.1) 


with  zero  flux  condition  and  which  is  similar  to  the  system  of  equations  in  (5.7)  without  external  noise 
term.  The  uniqueness  of  the  solutions  follows  from  the  uniqueness  of  this  system,  the  properties  of  linear 
and  bilinear  operators  and  by  using  the  Poincare  inequality  for  the  admissible  channel  domain. 

The  uniqueness  of  the  pressure  field  p(x,  t,  uj)  follows  from  the  uniqueness  of  the  constructed  basic 
scalar  field  P(x,  t,^;)  in  Sect.  3  and  the  uniqueness  of  the  perturbed  pressure  g'(x,  t,a;).  □ 


Remark  9.2.  Theorem  9.1  can  be  extended  to  the  case,  where  the  original  problem  (2.2)-(2.6)  has  a 
multiplicative  Gaussian  noise  or  a  multiplicative  Levy  noise  as  external  forcing.  For  this  case  the  basic 
vector  field  can  be  constructed  in  the  same  way  as  discussed  in  this  paper.  The  method  of  proof  of  the 
existence  and  uniqueness  of  the  perturbed  vector  field,  under  the  suitable  assumptions  (growth  property 
and  Lipschitz’s  condition)  on  the  noise  coefficients,  can  be  obtained  from  [44]. 
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